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ABSTRACT 

We present the results of 16 years of monitoring stellar orbits around the massive black hole in 
center of the Milky Way using high resolution near-infrared techniques. This work refines our previous 
analysis mainly by greatly improving the definition of the coordinate system, which reaches a long- 
term astrometric accuracy of ~ 300 /xas, and by investigating in detail the individual systematic 
error contributions. The combination of a long time baseline and the excellent astrometric accuracy 
of adaptive optics data allow us to determine orbits of 28 stars, including the star S2, which has 
completed a full revolution since our monitoring began. Our main results are: all stellar orbits are fit 
extremely well by a single point mass potential to within the astrometric uncertainties, which are now 
« 6x better than in previous studies. The central object mass is (4.31 ± 0.06|stat ± 0.36|i{o) x IQ^Mq 
where the fractional statistical error of 1.5% is nearly independent from i?o and the main uncertainty 
is due to the uncertainty in Rq. Our current best estimate for the distance to the Galactic Center is 
i?o = 8.33 ± 0.35 kpc. The dominant errors in this value is systematic. The mass scales with distance 
as (3.95 ± 0.06) x lO®(i?o/8kpc)^'^^M0. The orientations of orbital angular momenta for stars in the 
central arcsecond are random. We identify six of the stars with orbital solutions as late type stars, and 
six early-type stars as members of the clockwise rotating disk system, as was previously proposed. We 
constrain the extended dark mass enclosed between the pericenter and apocenter of S2 at less than 
0.066, at the 99% confidence level, of the mass of Sgr A*. This is two orders of magnitudes larger 
than what one would expect from other theoretical and observational estimates. 
Subject headings: blackhole physics — astrometry — Galaxy: center — infrared: stars 



1. INTRODUCTION 

Observations of Keplerian stellar orbits in the Galac- 
tic Center (GC) that revolve in the gravitational poten- 
tial created by a highly concen trated mass of ro ughly 
4 X 10^ Mp. JSchodel et al. 2002; Eis enhauer et al.ll2005l : 
iGhez et all [2003. 2005 ) currently constitute the best 
proof for the existence of an astrophysical massive black 
hole. In this experiment the stars in the innermost arc- 
second (the so-called S-stars) of our galaxy are used as 
test particles to probe the potential in which they move. 
Unlike gas, the motion of stars is determined solely by 
gravitational forces. Since the beginning of the obser- 
vations in 1992 one of the stars, called S2, has now 
completed one fu ll orbit. Its orbit (jSchodel et al.l[200l : 
IGhez et al. l'2003l has a period of 15 years. Since 2002 the 
number of reasonably well-determined orbits has grown 
from one to 28; in total we currently monitor 109 stars, 
see Figure [if . 

Due to the high interstellar extinction of w 30 magni- 
tudes in the optical towards the GC the measurements 
have to be performed in the near infrared (NIR), where 
the extinction amounts to only « 3 magnitudes. The 
first positions of S-stars were obtained in 1992 by Speckle 
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imaging at ESO's NTT in La Silla, a 4m-telescope, and in 
1995 a t the Keck telescop e, a lOm-telescope. Sin ce 1999 
Keck: IGhez et all (12001 ) and 2002 rVLT: iSchodel et all 
2002f )) the combination of 8m/10m-class telescopes and 
adaptive optics (AO) has been routinely used for deep 
(Hw 19) diffraction hmited (FWHM 40-100 mas) imag- 
ing and spectroscopy. 

The GC is a uniquely accessible laboratory for explor- 
ing the interactions between a massive black hole (MBH) 
and its stellar environment. By tracking the orbits of 
stars close to the MBH one can gather information on 
the gravitational potential in which they move. Of prime 
interest is the value of i?o, the distance to the GC, as it 
is one of the fundamental quantities in models for our 
Galaxy. Equally interesting is the nature of the mass 
responsible for the strong gravitational forces observed. 
While the measured mass makes a compelling case for a 
MBH, the exact form of the potential encodes answers to 
many interesting questions. Clearly, testing general rela- 
tivity for such a heavy object is among the goals; the first 
step would be to detect the Schwarzschild precession of 
the pericenters of some orbits. A measurable deviation 
from a point mass potential would give access to a possi- 
ble cluster of dark objects around the MBH, testing many 
theoretical ideas, such as mass segregation or the concept 
of a loss-cone. Another focus of interest are the proper- 
ties of the stellar orbits. The distributions of the orbital 
elements may have conserved valuable information about 
the formation scenario of the respective stars. This ad- 
dresses for example the so-called ' paradox of youth' for 
the stars in the central arcsecond (|Ghez et al.ll2003f) or 
the puzzling existence of a large number of 0-stars and 
Wolf-Rayet stars in the GC (Paumard ct al. 2006). 
This paper is the continuation of our long-term work on 
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Fig. 1. — Finding chart of the S-star cluster. This figure is based on a natural guide star adaptive optics image obtained as part of 
this study, using NACO at UT4 (Yepun) of the VLT on July 20, 2007 in the H-band. The original image with a FWHM of !a 75 mas 
was deconvolved with the Lucy-Richardson algorithm and beam restored with a Gaussian beam with FWHM= 2 pix=26.5 mas. Stars as 
faint as nifj = 19.2 (corresponding roughly to mx = 17.7) are detected at the 5cr level. Only stars that are unambiguously identified in 
several images have designated names, ranging from SI to S112. Blue labels indicate early-type stars, red labels latc-typc stars. Stars with 
unknown spectral type are labelled in black. At the position of Sgr A* some light is seen, which could be either due to Sgr A* itself or due 
to a faint, so far unrecognized star being confused with Sgr A*. 



stellar motions in the vicinity of Sgr A*. We reanalyzed 
all data available to our team from 16 years. The basic 
steps of the analysis are: 

• Obtain high quality, astrometrically unbiased maps 
of the S-stars. Obtain high quality spectra for these 
stars. 

• Extract pixel positions from the maps and radial 
velocities from the spectra. 

• Transform the pixel positions to a common astro- 
metric coordinate system; transform the radial ve- 
locities to the local standard of rest (LSR). For the 
astrometric data several steps are needed: 

— Relate the fainter S-stars positions to those of 
the brighter S-stars (Speckle data only). 

— Relate the S-stars positions to a set of selected 
reference stars. 



— Relate the reference stars to a set of SiO 
maser stars, of which the positions relative to 
Sgr A* are known with g ood accuracy fro m 
radio (VLA) observations ()Reid et al.ll2007f) . 

• Fit the data with a model for the potential and 
gather in that way orbital parameters as well as 
information about the potential. 

We organize this paper according to these steps. 

2. DATA BASE 

The present work relies on data obtained over many 
years with different instruments. In this section we 
briefly describe the different data sets. 

2.1. SHARP 

The first high resolution imaging data of the GC region 
were obta ined in 1992 with the SHARP camera built at 
MPE (Hofmann et"aL|[l99^ lEckart et al.|[l994 ). SHARP 
was used by MPE scientists until 2002 at ESO's 3.5m 
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NTT in Chile. The data led t o the detection of hig h 
proper motions close to Sgr A* (|Eckart fc Genze]|ll996f ). 
The camera was operating in speckle mode with exposure 
times of 0.3 s, 0.5 s and 1.0 s, which was the optimum 
compromise between sufficient signal to noise and fast 
sampling of the atmosphe r ic tur bulence. The data are 
described in lSchodel etall (|2003f ). a summary is given m 
table[Tl \ye use d the simple shift-and-add (ssa) technique 
(jChistou II1991D in order to obtain deep diffraction lim- 
ited images fr om the raw frames. Compared to our previ- 
ous analysis fjSchodel et al.ll20d3tlEisenhauer et al.ll2005D 
we did not base the astrometry on images combined from 
multiple pointing positions. Due to the camera's image 
distortions one should not trust the larger scale astrome- 
try of such multi-pointing images since coadding different 
pointings in the presence of static image distortions will 
lead to discontinuities in the effective distortion map of 
the combined image. These in turn would not be de- 
scribed well by the polynomial relations we use to map 
pixel positions onto the astro metric coordinate system, 
resulting in astrometric biases ()Reid et al.ll2003f ). Single 
pointing combinations are astrometrically unbiased - al- 
though not necessarily linear. They do not show disconti- 
nuities and represent smooth coordinate grids. Hence, we 
co-added the frames pointingwise, yielding typically four 
coadded maps per observing epoch ([Ott 2002). In order 
to get deeper in the central arcsecond, we also co-added 
all frames into one single map per epoch. Of the latter 
map the astrometry however can only be trusted for a re- 
gion as small as the central arcsecond (the region which 
was present in all pointing positions), leading to an addi- 
tional step of cross-calibration between the deep map and 
the four single-pointing maps per epoch. The ssa maps 
had a typical Strehl ratio of 30%. We further cleaned 
them using careful deconvolution and beam restoration, 
following the strategy outlined in Schodcl ct al. (2003). 
In order to assess the errors introduced by the decon- 
volution we used two different deconvolutio n methods: 
The Lucy-Richardson algorithm (iLucv II1974D and an it- 
erative blind deconvolution process l Jefferies fc Christoul 
[1993) , yielding two different (although not independent) 
maps. 

2.2. NACO 

The first AO imaging data available to us of the GC re- 
gion was obtained in 2002 with the Naos-Conica (NACO) 
system mounted at the fourth unit telescope Y epun of 
the VLT (jLenzen et al.iri998HRousset et al.lll998[ ). Com- 
pared to the SHARP data the NACO data are superior 
due to the larger telescope aperture (8.0 m versus 3.5 m) 
and the higher Strehl ratios (typically 40% for NACO) of 
the AO which is equipped with an IR wavefront sensor, 
allowing the use of the nearby K=6.5 mag star IRS7 as 
AO guide star. Furthermore, the sampling is increased 
compared to the Speckle data. For NACO we have typ- 
ically ten epochs per year, compared to one per year for 
SHARP. We obtained images both in the 27mas/pix and 
the 13mas/pix image scales. 

• In order to measure the positions of the SiO maser 
stars IRS9, IRSIOEE, IRS12N, IRS15NE, IRS17, 
IRS19NW, IRS28 and SiO-15 (|Reid et all 12001 . 
we used the 27mas/pix image scale both in H- 
and K-band in all years since 2002. The data 



are descri bed in iReid et al] ()2007f ): iTrippe et al.l 
(jin prep. I ) and summarized in table [21 The typ- 
ical single detector integration time was two sec- 
onds, such that the bright IR sources present in 
the r K, 20" field covered did not get saturated. 
Mostly, we used a dither pattern of four positions 
that guaranteed that the central ten arcseconds 
are imaged in each pointing position. The num- 
ber of useful maser positions per image varied be- 
tween six and eight. IRS19NW was not in the im- 
ages in 2002, 2003 and 2006; SiO- 15 was not cov- 
ered in 2003. Due to their brightness IRS17 and 
IRS9 were in the non-linear regime of detector in 
the observations from June 12, 2004 and thus ex- 
cluded for that epoch. Since the NACO camera 
when operated in the 27 mas/pix mode exhibits no- 
table geometric image distortions we constructed 
de-distorted mosaics from the individual images by 
applying a distortion correction, involving rebin- 
ning of the measured flux distribution to a new 
pixel grid. The proced ure is described in detail in 
iTrippe et al.l (|in prep.[ ) and relies on comparing dis- 
tances between stars present in the different point- 
ings. The distortion model used \sp = p'{l~f3p''^) 
with /3 « 3 X 10~^ where p and p' denote true and 
distorted pixel positions with respect to some ori- 
gin in the image that also is determined from the 
data. See also fig. [H We did not apply deconvolu- 
tion techniques on these images. 

• The positions of the S-stars were determined 
mostly from images obtained in the 13 mas/pix im- 
age scale. (Only when no image in the 13 mas/pix 
scale was available sufficiently close in time, we 
used also images obtained in the 27 mas/pix scale.) 
A typical data set contains two hours of data. The 
single detector integration time was mostly around 
15 seconds, and the field of view was moved after a 
few integrations successively to four positions such 
that the central four arcseconds are present in all 
frames. The data are summarized in table [1] and 
a complete list of the data sets used is given in 
the table in appendix [B] The reduction followed 
the usual steps of sky subtraction and flat-fielding. 
Manually selected high quality frames were com- 
bined to a single ssa map per epoch since the opti- 
cal distortions are s mall enough to be negl ected in 
the 13 mas/pix scale (jTrippe et aLllin prep.l ) for the 
frame combination. A distortion model of the same 
type as for the 27 mas/pix scale images was con- 
structed for each epoch; however the best-fitting 
model parameters varied more than expected be- 
tween the different epochs. We concluded that we 
were not able to solve for the distortion parame- 
ters with our observations. Hence, we did not ap- 
ply distortion models to the 13 mas/pix data but 
used higher order transformations when relating 
pixel positions to astrometric positions (see fig. [5]). 
In order to separate sources we moderately decon- 
volved the central five arcseconds of these maps 
with the Lucy-Richardson algorithm. The latter 
used a point spread function constructed from the 
map itself obtained b y applying the starfinder code 
([Diolaiti et al.|[2000[ ). In order to estimate the de- 
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TABLE 1 

Summary of the yearly number of epochs for which we 
obtained s-star images and the yearly mean number of 
s-star positions determined per epoch. 



Year Instrument # epochs # — 



1992 


SHARP 


1 


33 


1994 


SHARP 


1 


41 


1995 


SHARP 


1 


38 


1996 


SHARP 


2 


38.5 


1997 


SHARP 


1 


38 


1998 


SHARP 


1 


33 


1999 


SHARP 


1 


39 


2000 


SHARP 


1 


38 


2000 


GEMINI 


1 


31 


2001 


SHARP 


1 


39 


2002 


SHARP 


1 


21 


2002 


NACO 


8 


81.9 


2003 


NACO 


12 


83.8 


2004 


NACO 


11 


75.3 


2005 


NACO 


5 


86 


2006 


NACO 


13 


72.8 


2007 


NACO 


13 


101.2 


2008 


NACO 


5 


104.8 



TABLE 2 

Summary of the number of available maser star mosaic 

IMAGES, number OF MASER STARS PRESENT IN EACH FRAME AND 
THE RESPECTIVE FWHM OF THE POINT SPREAD FUNCTION IN THE 
IMAGES. 



Date # mosaics FWHM [mas] # maser stars 



May 2002 1 70 7 

May 2003 1 74 6 

June 2004 3 70, 81, 86 6, 6, 8 

May 2005 1 88 8 

April 2006 2 100, 100 7, 7 

March 2007 2 80, 80 8, 8 

March 2008 1 84 8 



convolution error we divided each 13mas/pix data 
set into two and obtained two coadded maps, each 
with half of the integration. Both maps were then 
deconvolved the same way as the full coadd. 

2.3. SINFONI 

Spectroscopy enables one to determine radial velocities 
of stars if the positions of known atomic or molecular 
lines can be measured in the stellar spectra. The GC is 
best exploited with integral field spectroscopy as one is 
interested in the radial velocities of all stars for which 
one can hope to determine orbits, i.e. all stars in the 
central arcsecond. In the NIR the K-band (2.0 — 2.4 /im) 
is best suited since it contains the hydrogen line Bracket- 
7 at 2.16612 /im. This line is present in absorption for 
B typ e stars, the most comm on spectral type for the S- 
stars (lEisenhauer et al1l2005f) . For late- type stars the 
CO band heads at 2.2935 /im, 2.3227 /im, 2.3535 /im and 
2.3829 /im are also covered by the K-band. 

Since July 2004 we regularly monitored the GC 
with the AO assisted field spectrometer SINFONI 
(jEisenhauer et all l2003t iBonnet et all l2004[ ). The in- 
strument is mounted at the Cassegrain focus of ESO's 
UT-4 (Yepun) and offers several operation modes con- 
cerning pixel scale and wavelength coverage. For the 



TABLE 3 

Summary of SINFONI data used for this work. The 

EXPOSURE TIME IS THE EFFECTIVE SHUTTER-OPEN TIME ON S2, 

FOR OTHER STARS THE ACTUAL EXPOSURE TIME MIGHT BE 
DIFFERENT SINCE THE OBSERVATIONS WERE MOSAICING AROUND 

Sgr a*. The FWHM was determined from a median image 

OF THE RESPECTIVE CUBE ON THE UNCONFUSED STAR S8. 



Date 


Band 


texp on S2 


FWHM 


# S-sts 






[min] 


[mas] 


with velo 


July 14 2004 


H-l-k 


40 


79 


7 


July 17 2004 


K 


110 


93 


25 


August 18/19 2004 


K 


80 


88 


23 


February 26 2005 


K 


20 


108 


4 


March 18 2005 


K 


10 


150 


4 


March 19 2005 


K 


40 


69 


16 


June 15 2005 


K 


200 


113 


8 


June 17 2005 


K 


440 


88 


25 


Aug 28 - Sep 5 2005 


K 


10 


>250 


5 


October 2-6 2005 


H+K 


120 


74 


22 


March 16 2006 


H+K 


110 


76 


27 


April 21 2006 


H+K 


10 


100 


6 


August 16/17 2006 


H-l-K 


100 


88 


18 


March 26 2007 


H-l-K 


20 


86 


10 


July 18-23 2007 


H+K 


133 


78 


15 


September 3/4 2007 


H+K 


70 


81 


15 


April 4-9 2008 


H+K 


200 


65 


40 


June 4 2008 


H+K 


10 


84 


3 



GC we operated SINFONI mostly in the AO scale, map- 
ping 0.8" xO. 8" onto 64x32 spatial pixels. We used the 
K-band grating and the combined H-l-K grating of SIN- 
FONI, with spectral resolutions of 4000 and 1500 respec- 
tively. For most of the data sets, the single exposure 
time per frame was 10 minutes; a few data sets also used 
5 minute exposures. We chose various mosaicking pat- 
terns inside the central arcsecond for the different runs; 
mainly with the aim to have a good compromise between 
monitoring the activity of Sgr A* and building up inte- 
gration on the S-stars. For stars at somewhat larger radii 
(r > 1") where confusion is less severe we also used data 
originally obtained for other scientific programs in the 
lOOmas/pix scale offering a field of view of 3.2"x3.2". 

The SINFONI AO works in the optical. Since the GC 
region is heavily extincted, one has to use a guide star 
relatively far away from Sgr A*. It is located 10.8" East 
and 18.8" North of Sgr A* and has a magnitude of mn = 
14.65. As a result the performance of the AO strongly 
depends on the seeing conditions. Therefore the quality 
of our SINFONI data is variable over the data set. For a 
typical run, one can detect Bracket-7 absorption of early- 
type stars as faint as mx = 15.5 and the CO band heads 
of late-type stars up to mx = 16.0. A summary of our 
data is given in Tabled 

We applied the standard data reduction for SINFONI 
data, including detector calibrations (such as bad pixel 
corrections, flat-fielding and distortion corrections) and 
cube reconstruction. The wavelength scale was cali- 
brated by means of emission line lamps and finetuned on 
the atmospheric OH lines. The remaining uncertainty of 
the wavelength scale corresponds to typically < lOkm/s. 
We did not trust the SINFONI cubes for their astro- 
metric precision, they were used only for their spectral 
dimension. Nevertheless it is easy to identify stars in the 
cubes. 



Stellar orbits in the Galactic Center 



5 



2.4. Other 

Beyond the data sets described so far, we included a 
few more data points which we describe briefly in this 
section. 

• Positions from public Gemini data for 2000: 

In addition to our observations we included images 
from the Galactic Center Demonstration Science 
Data Set obtained in 2000 with the 8-m-telescope 
Gemini North on Mauna Kea, Hawaii, using the 
AO system Hokupa'a in combination with the NIR 
camera Quire. These images were processed by the 
Gemini team and released to be used freely. We 
treated this data in the same way as the SHARP 
data. 

• Published radial velocities of S2 in 2002: The 

first radial measu rements of S2 were obtained by 
iGhez et all (|200^. We included the two pubhshed 
radial velocities since they extend the sampled time 
range by one year and clearly contribute signifi- 
cantly in fixing the epoch of pericenter passage tp 
for S2. 

• Radial velocities from longsplit spectroscopy 
with NACO in 2003: We used NACO in its 
spectroscopic mode to measure the radial veloc- 

ity of S2 in 200 3. The data are described in 

lEisenhauer et all (|2003f) . 

• Radial velocities from integral field spec- 
troscopy with SPIFFI in 2003: SPIFFI is the 
integral field spectrometer inside SINFONI. We 
used it without AO in 2003 as guest instrument 
at ESO-VLT UT-4 (Yepun) under superb atmo- 
spheric conditions and obtained cubes from which 
radial velocities for 18 stars (namely SI, S2, S4, SB, 
SIO, S12, S17, S19, S25, S27, S30, S35, S65, S67, 
S72 S76, S83, S95 S96) The data are described in 
lEisenhauer et all (|2003h . 

3. ANALYSIS OF ASTROMETRIC DATA 

This section describes in some technical detail the as- 
trometric calibration of our data. The first step is to mea- 
sure the positions of stars on the astrometric maps. Next, 
these positions of stars on the detector have to be trans- 
formed into a common astrometric reference frame. This 
procedure ultimately relies on measurements of eight SiO 
maser stars of which positions can be determined both 
in the radio and in NIR images. However, a direct com- 
parison of the central arcsecond and the maser stars on 
one and the same image is impractical for two reasons: 
a) The exposure times necessary to obtain sufficiently 
deep images for the S-stars saturates the detector at the 
positions of the maser stars, b) The field of view of the 
13 mas/pix pixel scale is too small to show enough maser 
stars. Therefore we need to crosscalibrate the S-stars im- 
ages with the maser star images. This is done by a set of 
selected reference stars (fig. which are present both 
in the S-star images and the maser star images. For the 
SHARP data, even an additional step of cross-calibration 
is taken. We selected reference stars with 1" < r < 4" 
that are brighter than mx ~ 14.5 and apparently uncon- 
fused, yielding a sample size of 91 stars. 



3.1. Extraction of pixel positions 

All pixel positions were obtained by two-dimensional 
Gaussian fits in the images. The fits yielded both the 
positions and estimates for the statistical error of the 
positions (section I3.4.3P . For each epoch for which we 
have useful S-stars data we extracted pixel positions for 
the S-stars and for the reference stars. 

3.1.1. SHARP 

Only star images that are not visually distorted (e.g. 
due to a confusion event) were used from the SHARP 
data. 

• Reference stars: We obtained the reference stars' 
positions from the four single-pointing maps from 
each epoch. Due to the limited field of view, in 
each frame only a subset of the reference stars is 
present. 

• Brighter S-stars: For the brighter S-stars (e.g. 
S2, SI, S8, SIO, S30, S35) typically aU four differ- 
ent pointing positions could be used. The astro- 
metric position of each star was determined from 
the corresponding four pixel positions using the as- 
trometric average position (see section I3.4.3P . 

• Fainter S-stars: In order to detect faint S-stars 
we used the fifth coadded map which can be trusted 
astrometrically only for the innermost arcsecond. 
The limiting magnitude for a non-confused source 
was typically tok ~ 15.8. We determined the pixel 
positions of the weaker S-stars as well as the ones 
of the brighter S-stars. The latter served as refer- 
ence for relating the fainter stars to the astrometric 
coordinate system (see section |3. 4. 3|) . 

Since we had two different deconvolutions at hand, we 
extracted pixel positions from both sets of images. Thus, 
up to eight (= two deconvolutions x four pointings) pixel 
positions were obtained per star and epoch. 

3.1.2. NACO 

For the NACO data, we used both the 27 mas/pix data 
and the 13 mas/pix data. 

• SiO maser stars: Positions for the SiO maser 
stars were obtained by Gaussian fits to the stars' 
images in the 27 mas/pix mosaics. The SiO maser 
stars were unconfused in all mosaics. 

• Reference stars: The positions of the reference 
stars were measured both on the 27 mas/pix mo- 
saics as well as on the 13 mas/pix maps (Table [Ij, 
since they serve as cross-calibration between the 
two sets. They were selected to be unconfused, 
thus essentially it was possible to use all reference 
stars visible on any given frame. 

• S-stars: For isolated S-stars, the positions were 
obtained from a simple Gaussian fit to the manu- 
ally identified stars in the maps. Due to the higher 
sampling rate with NACO confusion events can 
be tracked much better in the AO data than in 
the SHARP data. Therefore it was reasonable to 
also measure positions when stars are partly over- 
lapping. In such simultaneous, multiple 
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Gaussian fit to the individual peaks was used, re- 
sulting of course in larger statistical uncertainties 
of the obtained positions. 
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Fig. 2. — The open symbols mark the sample of 91 reference 
stars which are used to define the astrometric frame for the S-stars. 
The underlying image was obtained on April 3, 2007 in H-band, 
deconvolved and beam-restored with a beam of 2 pix. North is up, 
East is left. The field is 9.3" x 9.3". 



3.2. Relating the reference stars to the SiO maser stars 

The goal of this step is to obtain linear models for 
the motions of the reference stars, i.e. to express their 
velocities and positions with linear functions x{t),y{t) 
in terms of astrometric coordinates. These models then 
define a common reference frame that is calibrated in 
position and velocity such that radio Sgr A* should be at 
rest at the origin of the system. Such a coordinate system 
allows one to test if the center of mass obtained from 
orbital fitting coincides with the compact radio source^. 

In the following we present two ways to obtain the de- 
sired calibration. They differ in the way in which the 
positions and velocities for the reference stars are deter- 
mined: either all maser star images are tied to the re- 
spective maser positions (multi-epoch cross calibration); 
or only one maser star image is used to tie to the ra- 
dio maser positions, and the other maser star images are 
matched to that by an additional step of cross calibration 
that only involves infrared data (single epoch cross cali- 
bration). It turns out that both ways have their specific 
advantages and disadvantages in terms of position and 
velocity calibration of the resulting coordinate systems. 
We finally constructed a third coordinate system com- 
bining the advantages and rejecting the disadvantages. 

3.2.1. Multi-epoch cross calibration with all maser star 

images 



Systematic problems of the coordinate system could be ab- 
sorbed into the orbital fitting by allowing the center of mass to 
have an offset from 0/0 and a non-zero velocity, at the cost of not 
being able to test the coincidence of the center of mass with radio 
Sgr A*. 



Using the results from lReid et al.l (|2007D we calculated 
the expected radio maser positions for the given obser- 
vation epochs. The different maser images contained be- 
tween seven and nine SiO masers of which we used six 
to eight since we excluded IRS7 due to its brightness of 
niK ~ 6.5. By allowing for a linear transformation of 
type X = xq + M.p between the astrometric positions x 
and the pixel positions p in the respective image we de- 
termined a transformation by which any detector posi- 
tion can be converted into astrometric coordinates. Note 
that the use of a linear transformation is justified since 
the IR images were distortion-corrected mosaics. The 
rms of the ID-residuals of the SiO masers (thus apply- 
ing the transformation to the SiO masers' pixel positions 
and comparing the result with the expected radio po- 
sitions) was 2.28 mas. Correspondingly we expect that 
from our 11 images a coordinate system can be defined 
to at most an ID-accuracy of 2.28/vTTrnas« 0.7mas if 
the measurement errors from the 11 images are uncor- 
related. The transformation was applied to the sample 
of reference stars in each image. We then fitted the re- 
sulting astrometric positions of the reference stars as a 
function of time with linear functions. From these lin- 
ear fits we obtained residuals, allowing obvious outliers 
to be identified and rejecting them. We excluded stars 
that had a residual different from the median residual by 
twice that value, if the deviation was larger than 2 mas. 
That excluded between and at most 10 of the 91 ref- 
erence stars for the various mosaics, the reason for these 
outliers being confusion that in some mosaics affects the 
fainter reference stars due to the varying image quality. 
After this moderate data cleaning we repeated the fit- 
ting. The mean rms of the ID-residuals per image had 
then a value of 1.45 mas. 

The next step of refinement was to compare all mea- 
sured positions in one mosaic with the positions expected 
from the fits, effectively checking how well a given im- 
age fits to the other 10 images. A visual inspection of 
maps of residual vectors showed that the residuals are 
not randomly distributed but unveiled some systematic 
shift and rotation for each image. Since each image is 
compared with 10 other images, any systematic prob- 
lem in the given image is most likely to come from that 
image and not from a combined effect of the others. In- 
deed, the interpretation of the observed systematic effect 
is straightforward, it means that each individual mosaic 
is not registered perfectly with respect to the sample av- 
erage, i.e. the transformation for the respective image 
is slightly wrong. This systematic error is naturally ex- 
plained by measurement errors of the positions of the 
SiO maser stars in the respective image. Such an er- 
ror translates into an error of the parameters of the lin- 
ear transformation used to tie the astrometric frame to 
the pixel positions in the mosaic and shows up as a sys- 
tematic effect in the residuals of the independent set of 
reference stars. Thus, we were able to determine bet- 
ter transformation parameters by adding to the original 
linear transformation the linear transformation that min- 
imizes the residuals of the reference star sample, yielding 
a corrected linear transformation. We applied it to the 
data and obtained the final linear motion models for the 
reference stars. The rms of the ID-residuals now was 
0.55 mas. This step changed the position of the origin 
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by (Aa, A6) = (-0.01, 0.05) mas and the velocity of the 
system by less than 4/xas/yr, these quantities being the 
mean differences of the respective quantities for the ref- 
erence stars before and after the refinement. Hence, the 
refinement effectively did not change the coordinate sys- 
tem calibration. We call the coordinate system so defined 
the 'maser system' in the following. 

The position of the origin of the maser system and its 
velocity are uncertain due to two effects: a) the non- 
zero errors of the SiO maser stars' radio positions and 
velocities and b) the IR positions of the SiO maser stars 
show some residuals to the best fitting linear motion, 
indicative of residual image distortions and of measure- 
ment errors in the pixel positions in the IR images. The 
propagation of the statistical errors into the definition 
of the coordinate system was addressed using a Monte- 
Carlo technique. We varied the input to the transfor- 
mations according to the measured errors and residuals. 
We created 10^ realizations of transformations, assuming 
a Gaussian distribution of the simulated values around 
the original values. The standard deviation of the posi- 
tions obtained for Sgr A* estimates the positional uncer- 
tainty of the maser system under the assumption of un- 
correlated measurement errors. We obtained (Aa, AS) = 
(0.46, 0.77) mas. Similarly, the standard deviation of the 
velocities obtained for Sgr A* estimates the uncertainty 
of the maser system's velocity under the same assump- 
tion. We obtained {Ava,Avs) = (0.29, 0.55) mas/yr. 

However, in our data the assumption that the errors 
from the 11 maser images are uncorrelated is not fulfilled. 
We rather observe a typical residual per SiO maser star 
for all epochs when comparing the transformed, mea- 
sured positions with the predicted radio positions. Pos- 
sible reasons are: first, the linear motion models obtained 
for the SiO masers could be inaccurate due to some un- 
known some unknown systematic problem of the radio 
positions. Secondly, the radio positions could not be ap- 
plicable to the IR positions, for instance if the maser 
emission would originate from far away of the stellar sur- 
face. Thirdly, the correlation could arise due to some 
unaccounted systematics in the infrared frames, such as 
uncorrected distortion. We fitted the residuals of each 
star with linear functions and obtained in that way es- 
timates for the mean position and mean velocity uncer- 
tainty for each star. Then we calculated the mean devi- 
ation (over the SiO maser stars which are < 15" away 
from Sgr A*) of these linear motion model parameters 
as estimates for the positional and velocity uncertainty 
of the maser system given the correlations in our data. 
With our initial transformation we obtained {Aa, A6) = 
(0.92 ± 0.42, 2.22 ± 0.43) mas and {Ava,Avs) = (0.41 ± 
0.24, 0.29 ± 0.24) mas/yr. After the refinement we 
got (Aa, AS) = (0.95 ± 0.73, 2.35 ± 0.58) mas and 
{Ava,Avs) = (0.38 ±0.41, 0.28 ± 0.33) mas/yr. Finally 
we conservatively adopt for the uncertainties of the maser 
system {Aa,AS) = (1.0, 2.5) mas and Ava — Avs = 
0.5 mas/yr. The positional uncertainty is considerably 
larger than what one would have obtained for uncorre- 
lated residuals. 

While the maser system is a direct crosscalibration of 
maser and reference stars, the resulting velocities of the 
reference stars are directly sensitive to errors in the ve- 
locities of the SiO maser stars, both in the radio data 
and the NIR 27mas/pix mosaics. 



3.2.2. Single-epoch cross calibration with one maser star 

image 

The sensitivity of the reference star velocities to er- 
rors in the SiO maser velocities can be avoided by an 
additional step of cross-calibration: We can measure the 
positions of the reference stars in all maser images with 
respect to a much larger sample of stars in these images. 
This cluster of stars is assumed to be non-rotating and 
not moving with respect to Sgr A*. The cluster is tied to 
the astrometric frame for just one epoch, as given by the 
radio positions of the SiO mase rs, which c a n be calcu- 
lated for the chosen epoch from iReid et al.l (|2007f ). For 
all other epochs it is assumed that the mean cluster is 
stationary in time. Hence, the velocity calibration relies 
on the statistical argument that for a sufhciently large 
sample of cluster stars the mean velocity of the cluster is 
expected to become very small. For a typical velocity of 
V for a cluster star and N stars, the error of this mean 
should be of order u/ViV. 

Effectively the few maser stars are only used once in 
this scheme. Any error in their radio positions, radio ve- 
locities or NIR detector positions will therefore translate 
into a positional offset, but not into a systematic velocity 
of the coordinate system. The latter is instead connected 
to the validity of the assumption that the cluster mean 
is stationary. In order to ensure the best estimate of the 
velocity calibration we adopted the following procedure: 

1. We selected the maser mosaic from 12 May 2005, 
which was chosen since it is of good quality and 
roughly corresponds to the middle of the range in 
time co vered with NACO. Buil ding upon the work 
done bv lTrippe et al.l (|in prep.! ) we selected an en- 
semble of stars in that mosaic of which the positions 
can be measured with high reliability: Take all 
stars that have a peak flux of more than 25 counts 
which at the given noise level of 1.9 counts selects 
high-significance stars. In a second step many stars 
get excluded again: All stars with more than 700 
counts (they could be saturated in other frames 
with longer single detector integration times) and 
all stars that have a potential source (peak with 5 
counts) within 10 pixels. Furthermore a Gaussian 
fit was required to yield a FWHM< 0.05 pix and 
the fitted position must coincide with the position 
obtained using DAOPHOT FIND. This yielded a 
sample of 433 stars. We determined the astromet- 
ric positions at the given epoch of the 433 stars by 
means of a linear transformation that was deter- 
mined from the eight maser stars. 

2. We then determined preliminary astrometric posi- 
tions for a much larger sample of 6037 stars in all 
11 mosaics, by tying their pixel positions at all 11 
epochs to the astrometric positions of the 433 stars 
at the reference epoch with a linear transformation. 
Note that not taking into account proper motions 
in that step makes the velocity calibration indepen- 
dent from the radio measurements. The error due 
to the omission of the proper motions is minimized 
by using 433 stars instead of few masers for the 
cross-calibration. 

3. We fitted linear motion models to all 6037 stars. 
After the fit we determined the residuals of all star 
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positions in all mosaics. By inspecting the residu- 
als of any mosaic as function of position, we were 
able to map the residual image distortions in the 
given mosaic. These residual image distortions can 
arise due to imperfect registration of the individ- 
ual exposures to the respective mosaic, or due to 
an error in the distortion correction applied to the 
individual frames. 

4. For each star in each mosaic we determined an es- 
timate of the residual image distortion by calculat- 
ing the mean of the residuals of the stars in the 
vicinity (r < 2") of the chosen star. The radius 
was chosen such that a suitable number of stars 
was present in the area from which the correction 
was determined and the area was sufficiently lo- 
cal. A value of r < 2" was a good compromise, 
yielding 30-50 stars typically. The estimate for 
the residual image distortion was then subtracted 
from the given star; the typical values applied were 
Ax = 0.66 ± 0.20 mas and Ay = 0.62 ± 0.24 mas. 

5. In a second fit we used the corrected astrometric 
positions in order to obtain updated linear motions 
models for the 6037 stars. 

6. Then we defined the final cluster: it consists of all 
stars which were present in all 11 mosaics, with 
radii between 2" and 15" and which are not known 
early-type stars. These criteria yielded a cluster 
sample size of 2147 stars. 

7. We determined the cluster mean velocity, yielding 
(—0.04, 0.00) mas/yr, and subtracted that value 
from all velocities of the 6037 stars. This then is 
the final, velocity-calibrated list of linear motions 
from which the reference star sample is extracted. 
The mean radius of the 2147 stars of the cluster 
sample is 9.89" , the root mean square (rms) speed 
of the stars in the sample is 157km/s « 4.15 mas/yr 
(for Rq = 8kpc). 

We call the coordinate system defined in this way the 
'cluster system'. Since we expect the mean of the clus- 
ter to show a net motion of order 157/\/2147km/s — 
3.4km/s = 0.09 mas/yr, we estimate the uncer- 
tainty of the velocity calibration to be of the same 
size. We checked this number more thoroughly by 
means of a Monte Carlo simulation: We divided 
the cluster into nine radial bins with boundaries 
[2, 5, 7, 8.5, 10, 11, 12, 13, 14, 15] mas (selected such that in 
each bin roughly the same number of stars is present). 
For each bin we have determined the rms velocity 
yielding 3.7, 3.3, 3.0, 2.9, 2.8, 2.7, 2.6, 2.6, 2.6mas/yr re- 
spectively. We then simulated clusters in proper motion 
space, for each bin the Gaussian width of the velocity dis- 
tribution was set to the respective rms velocity and the 
number of stars was matched to the real numbers in each 
bin. For each simulated cluster we were able to obtain 
in that way a mean velocity; simulating 10000 clusters 
allowed us then to estimate the uncertainty of the mean 
cluster velocity. We obtained 0.06 mas/yr; even a bit bet- 
ter than the simple estimate. Hence, if the assumption 
of isotropy is correct, the cluster system should allow for 
a better calibration of the reference star velocities than 



with the maser system. The assumption could be wrong, 
for example if a net streaming motion were present in 
the GC cluster. 

The statistical positional uncertainty of the origin of 
the cluster system was estimated by the same means 
as for the maser system. We obtained (Aa, AS) — 
(0.85, 1.51) mas. In addition to these uncertainties, the 
residuals of the SiO masers also need to be considered, 
for the epoch at hand the mean deviation is (Aa, A6) — 
(1.87, 3.12) mas. The uncertainties here are greater than 
the respective numbers for the maser system due to the 
fact that the position of Sgr A* in the cluster system 
is effectively measured only on one frame while in the 
maser system it is measured in several and the residuals 
are not fully correlated. 

3.2.3. The final, combined coordinate system 

The maser system has a smaller systematic error in its 
position calibration, while the cluster system is superior 
with respect to the velocity calibration. Hence, by com- 
bining the two we were able to construct a system that 
combines both advantages. The idea simply is to correct 
cither the velocity calibration of the maser system such 
that it agrees with the one from the cluster system or to 
correct the origin of the cluster system such that it co- 
incides with the origin of the maser system (taking into 
account that the systems refer to two different epochs). 
Note that this implicitly uses the fact that the second, re- 
fining transformation of the maser system did not change 
its calibration properties. 

We used the sample of reference stars to compare the 
two systems. The mean positional offset between the two 
lists of positions for the epoch of the cluster system was 

/-1.87\ , /0.04\ 

PCSys-PMSys= K 37 I ± ( 0.O4) ' ^ ' 

Here, 'Csys' denotes the cluster system, 'MSys' the maser 
system. The errors are the standard deviation of the 
sample of differences. We also calculated the differences 
of the reference star velocities, as given by the two lin- 
ear motion models obtained for each reference star. We 
obtained 

/-0.60\ ^ /0.08\ mas 

VcSys - VMSys - ^^q^q j ± j — , (2) 

where again the errors are the standard deviation of the 
sample of differences. 

This means that the two coordinate systems differ sig- 
nificantly in position and velocity calibration in a sys- 
tematic way. It should be noted that only the difference 
between the two coordinate systems is that well defined; 
for the question how well each of the coordinate systems 
relates to Sgr A*, the larger, systematic errors of sec- 
tions 13.2.11 and 13.2.21 need to be considered. It is exactly 
the fact the difference between the coordinate systems 
is well defined that allowed us to combine the two co- 
ordinate systems and to gain accuracy in the combined 
system that way. Also note that the size of the offsets 
occurring here are consistent with the combined uncer- 
tainties of the two coordinate system; much larger offsets 
would have meant that the coordinate systems would be 
inconsistent with each other. 

Finally we chose the method which corrects the cluster 
system by a positional offset. The positional difference 
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from equation ([T]) was subtracted from all positions of the 
cluster stars (and thus also from the reference stars that 
are a subset of the cluster). This combined coordinate 
system has the same prior as the cluster system, namely 
that the cluster is at rest with respect to Sgr A*. The 
linear motion models so obtained were then used for the 
further analysis. 

3.3. Relating the S-stars to the reference stars 

We constructed the transformation from pixel posi- 
tions on the detector to astrometric positions by means 
of the reference stars. For each given image, we calcu- 
lated the expected astrometric positions of the reference 
stars using the linear motions models as obtained in sec- 
tion l3.2.3l Given the pixel positions of the reference stars 
in the respective image, we related the two sets of posi- 
tions by means of a cubic transformation (20 parameters) 
of type 

a;sky ^Pa+Pix+ P2V + Pax'^ + Pixy + p^^y^ + 

P&x^ + Pix'^y + psxy'^ + pgy^ 
J/sky ^qo + qix + q2y + q^x^ + qixy + q^y^ + 

q&x^ + qix'^y + qsxy'^ + qgy^ . (3) 
Once the transformations are known, it is straight- 
forward to apply them to the pixel positions of the S- 
stars. 

The parameters pi , qi were found by demanding that 
the transformation should map the two lists of positions 
optimally in a sense. Since the problem is linear, the 
parameter set can be found with a pseudo-inverse matrix 
(we always used at least 50 stars, thus 100 coordinates, 
for 20 parameters). The procedure also allows for an 
outlier rejection. For this purpose we applied the trans- 
formation to reference stars themselves and calculated 
the residuals to the expected astrometric positions. By 
only keeping reference star positions which are not more 
off than 15 mas from the expected position we cleaned 
our sample. This excluded in total 19 of the 7189 refer- 
ence star positions. For the cleaned set we redetermined 
the linear motion model for each star under the side con- 
dition that the refinement would not change the mean 
position or the mean velocity of the sample of reference 
stars, thus avoiding a change of the origin of the coor- 
dinate system and a change of its velocity. Compared 
to previous work the number of reference stars used is 
roughly a factor eight larger. This reduced the statisti- 
cal uncertainty of this calibration step to a very small 
leveF. 

For the SHARP data we had to use some additional 
steps for relating the S-star positions to the reference 
stars, since for a given epoch we used two deconvolu- 
tions for which we had four single-pointing frames and 
one combined map respectively. We used the pixel posi- 
tions of the reference stars in the two times four single- 
pointing images together with the predicted astrometric 
positions of the reference stars to set up eight transfor- 
mations of the type given in equation [3l Not all refer- 
ence stars are present in all pointings, but in all cases 

Actually some of the reference stars relatively close to Sgr A* 
were also considered as S-stars for which we tried to determine 
orbits. Indeed, four of those stars showed significant accelerations. 
However, we did not exclude them from the sample of reference 
stars. Therefore, an additional, obvious step of refinement would 
be to allow for quadratic motion models for the reference stars. 



their number exceeded 50, such that the transformation 
parameters were well determined. With these transfor- 
mations we calculated the astrometric positions of the 
brighter S-stars detected in the eight frames and used 
the average astrometric position in the end. The stan- 
dard deviation of the eight astrometric positions was in- 
cluded in the error estimate. For the fainter S-stars we 
used the coadded maps. For the two coadded maps (two 
deconvolutions) per epoch we set up two times four full 
first order transformations relating pixel positions of the 
brighter S-stars in each coadded map to the respective 
pixel positions in the four single-pointing frames. With 
these transformations we determined the pixel positions 
of the fainter S-stars which they would have had in the 
single-pointing frames. These fictitious pixel positions 
were then transformed with the cubic transformation of 
the respective single-pointing frame into astrometric po- 
sitions. The average of the latter was used in the end, the 
standard deviation was included in the error estimate. 

3.4. Estimation of astrometric errors 

The goal of this section is to understand the errors of 
the astrometric data. This includes both statistical and 
systematic error terms. The statistical error is due to the 
uncertainty of the measured pixel positions. Among the 
systematic error terms are the influence of the coordinate 
system, residual image distortions, transformation errors 
and unrecognized confusion. 

3.4.1. Offset and velocity of the coordinate system 

The accuracy in 2D-position {Ax, Ay) and 2D- velocity 
{Avx, Avy) of the combined coordinate system is given 
by the numbers in sections 13.2.11 and 13.2.21 In the third 
dimension, we don't use any priors for Az, since we wish 
to determine Rq from our data. 

For Avz we use the prior that Sgr A* is not mov- 
ing radially, based both on theoretical arguments and 
on radio and NIR measurements. Even if Sgr A* is dy- 
namically relaxed in the stellar cluster surrounding it, 
some random Brownian motion due to t he i nteraction 
with the surrounding stars is expected. [Merritt ct ahl 
(|2007l) calculated this number and concluded that the 
motion should be » 0.2km/s. This is consistent with 
the flndings of lReid fc Brunthaled ()2004[ ) who show that 
Sgr A* has a proper motion of = 18± 7km/s in galac- 
tic longitude and Vb — —0.4 ± 0.9km/s in galactic lat- 
itude (assuming Rq = Skpc). The significance of the 
fact that ^ is disputed, and furthermore it is not 
clear, whether it is truly due to a peculiar motion of 
Sgr A* or due to a difference between the global and local 
measure s of t he angular r otation rate of the Milky Way 
()Reid fc Brunt haleill2004f l. Clearly, the motion of Sgr A* 
perpendicular to the galactic plane is very small as ex- 
pected. In the third dimension, the velocity of Sgr A* 
can only be determined indirectly by radial velocity mea- 
surements of the stellar cl uster surrounding it. Using a 
sample of 85 late-ytpe stars lFiger et al.l ()2003f ) found that 
the mean radial velocity of the cluster is consistent with 
0: Vz ^ -10 ± llkm/s. iTrippe et all j in prep.f ) used 
a larger sample of 664 late-type stars and found con- 
sistently Vz — 4.6 ± 4.0km/s. Compared to that, the 
uncertainty AU « 0.5km/s in the defin ition of the lo- 
cal st andard of rest is much smaller (Dehnen fc Binnevl 
Il998f) . We conclude that all measurements are consistent 
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with Sgr A* being at rest at the dynamical center of the 
Milky Way and we assume a prior oi Vz = ± 5 km/s for 
our coordinate system. 

Summarizing, our combined coordinate system should 
be accurate to the numbers listed here, of which finally 
used the conservatively rounded values.: 

Aa; = 0.95 mas « 1.0 mas 

Ay = 2.35 mas w 2.5 mas 
Avx — 0.06 mas/yr « 0.1 mas/yr 
Avy = 0.06 mas/yr « 0.1 mas/yr 
At)^ = 5km/s. (4) 

3.4.2. Rotation and pumping of the coordinate system 

Potentially, there are two more degrees of freedom, 
which could affect the reliability of the chosen coordi- 
nate system, namely rotation and pumping. An artifi- 
cial rotation can be introduced if the selected stars by 
chance preferentially move on tangential tracks with a 
preferred sense of rotation. Similar, artificial pumping 
can occur: suppose that by chance all selected stars move 
on perfect radial trajectories and that stars further out 
move faster than stars closer to Sgr A*. Such a pattern, 
which would be somewhat similar to the Hubble flow of 
galaxies, would yield under the set of transformations 
a time-dependent plate scale and otherwise stationary 
stars. Both effects can affect the selection of the refer- 
ence star sample and (less important) the selection of 
cluster stars. 

The chosen coordinate system relies on the assumption 
that the cluster does not show any net motion (see sec- 
tion l3.2.2|) . net rotation or net pumping. The selection of 
a finite number of cluster stars however limits the accu- 
racy with which these conditions can be satisfied. Given 
2147 stars with a RMS velocity of « 157km/s and a typ- 
ical distance of 10" we expect that any selection leads to 
a pumping or rotation effect of the order of 9/ias/yr/". 

Due to the errors in the SiO maser positions, the maser 
system can show artificial pumping or rotation. Similar 
to what was done in sections f3.2.1l and l3.2.2l we simulated 
in a Monte Carlo fashion the error propagation. From 
10^ realizations of the transformations, assuming the ob- 
served errors of the SiO maser positions in the NIR and 
radio, we created perturbed sets of reference stars. The 
standard deviation of the pumping and rotation motion 
{Vr/r and Vt/r respectively) over these sets then estimate 
the stability of the maser system. We obtained 

Vr/r\MSys = ^7 ^J,a.s/yr/" , 

i't/?'|MSys = 33/xas/yr/" . (5) 

The cluster system (and therefore also the combined 
system) can be checked against the maser system. By 
calculating the difference in velocity for each reference 
star and subtracting from those the difference of the two 
coordinate system velocities we obtained a vector field of 
residual velocities, which is well described by: 

Vr/r\cSy - Wr/T-lMSy = (32 ± 2|stat ± 9|sys) Aias/yr/" 

vt/rlcsy - vt/r\MSy = (6 ± 2|stat ± 9|sys) //as/yr/" (6) 

The combined size of the effects from equations [5] and |6] 
estimate the error made when using the assumption that 
the combined coordinate system is non-rotating and non- 
pumping. At 1" these effects can sum up over 15 years 



to at most 0.7 mas, while for the maximum projected 
distance of S2 (« 0.2") the resulting positional errors 
are even a factor 5 smaller. We therefore neglected these 
effects in the following. 

3.4.3. Statistical errors of the pixel positions 

This paragraph deals with the uncertainties of the stel- 
lar positions on a given image; the unit of this error term 
as measured is therefore pixels. The error which is most 
easily accessible is the formal fit error of the Gaussian fit 
to a source. However, in deconvolved and beam-restored 
images it might be a bad estimator for the positional un- 
certainties. Therefore we compared additionally different 
deconvolutions of the same image for each epoch in order 
to get a more robust estimate. 
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Fig. 3. — The statistical errors of the pixel positions for the 
NACO K-band data as a function of arbitrary detector units of 
flux. The thin lines show the respective error model for each epoch; 
the thick dashed line is the mean for the data. The mean has a 
floor at 99 fias, the median (not shown) at 84 fias. 



For the SHARP data we used up to eight (= two decon- 
volutions X four pointings) pixel positions. The standard 
deviation of the astrometric positions was included in the 
error estimate for the statistical position error. For stars 
which were present only in one frame, the typical error 
of the epoch was used instead. 

For NACO we split up each data set into two parts and 
deconvolved both co-added images with the same point 
spread function as the co-added image of the complete 
data set (see Section [521) ■ We determined the pixel posi- 
tions of the reference and S-stars in the two deconvolved 
frames and applied a pure shift between the two lists 
of pixel positions such that the average pixel position 
is the same for both. The remaining difference between 
respective positions of one star estimates the statistical 
uncertainty for that star. The error estimates obtained 
this way were a strong function of the stellar brightness. 
Therefore we described the error estimates as a function 
of flux for each epoch (see Figure [3]) using a simple em- 
pirical model of the form ax~"' + b. The mean floor b over 
all data sets is 99/ias, while for lower fluxes the error in- 
creases up to 2 mas. We used the empirical description 
of each image to assign an error to all stellar positions 
obtained from that frame. Finally we checked whether 
the formal fit error of the positions was greater than the 
estimate from the empirical error model. In such a case 
we used the formal fit error instead. Figure H] shows 
the final distribution of statistical errors for the NACO 
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data. It is effectively the mean error model folded with 
the brightness distribution of the S-stars. 
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Fig. 4. — The measured distribution of the statistical errors of 
the pixel positions for the NACO data. The characteristic statis- 
tical error (defined as the peak of the distribution) is 108 /ias, the 
systematic error terms have to be added to this to come to a fair 
estimate of the true uncertainty. 
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Fig. 5. — Determination of residual image distortions for the 
NACO H-band data from September 8, 2007, 13mas/pix. The his- 
tograms show the differences of detector distances for a set of bona 
fide stars as measured in the four pointing positions with a dither 
offset of 7" . Left: Using the raw frames. Middle: After application 
of a distortion model, Right: After transforming the raw positions 
with a cubic transformation onto a common grid. The correspond- 
ing ID coordinate errors are determined from Gaussian fits to the 
distributions and are quoted at the top of each panel. 



For the SHARP data we obtained a broad distribution 
of the statistical pixel position errors with no clear max- 
imum and a tail to 2 mas. The median error is 360 /ias, 
the mean error 760 /las in the SHARP data. 

3.4.4. Residual image distortions 

A main source of error at the sub-milliarsecond level is 
image distortions. We estimated this error term by com- 
paring distances of stars in different pointing positions 
with a dither offset of 7" (see Figure [5]). If we had used 
only the raw positions and linear transformations, the re- 
sulting mean ID position error would be as large as 1 mas 
for the 13 mas/pix NACO data. By applying a distortion 
model (see section [2?2|) plus a linear transformation this 
error can be reduced to 600 /xas. Allowing for a cubic 
transformation onto a common grid yields an error of 
240 /ias only. This justifies our choice to use a high order 
transformation rather than to de-distort the 13 mas/pix 
NACO images. The numbers obtained in this way are 
actually the combined error of the statistical and trans- 
formation uncertainties with the residual image distor- 
tions. Subtracting the former we conclude that residual 
image distortions have a contribution of 210 /ias to the 
error budget of each individual astrometric data point. 
We thus added this value in squares to all other error 
terms, effectively acting as a lower bound for the astro- 
metric errors. 

We applied the same analysis to the 27 mas/pix NACO 
data which had a dither offset of 14" (see Figure [6]). The 
raw differences showed a skewed distribution, indicating 
the presence of image distortions. The rms of this dis- 
tribution is 2.1 mas. After applying the distortion model 
the typical residual error is reduced to 1.3 mas and the 
distribution is a nice Gaussian. Interestingly, mapping 
the positions with a cubic transformation onto each other 
does less well here. The distribution becomes less skewed, 
however it is still non-Gaussian and the rms is 1.6 mas. 
This justifies a-posteriori the use of the distortion cor- 
rection for the 27 mas/pix NACO data when determin- 
ing the motion models for the reference stars. For the 
SHARP data we obtained a characteristic error due to 
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Fig. 6. — Determination of residual image distortions for the 
NACO K-band data from March 16, 2007, 27 mas/pix. The his- 
tograms show the differences of detector distances for a set of bona 
fide stars as measured in the four pointing positions with a dither 
offset of 14". Left: Using the raw frames. Middle: After appli- 
cation of a distortion model. Right: After transforming the raw 
positions with a cubic transformation onto a common grid. The 
corresponding ID coordinate errors are quoted at the top of each 
panel, in the middle panel the value is determined from a Gaussian 
fit, for the other two the rms is quoted due to the non-Gaussianity 
of the distributions. 



residual image distortions of 0.8 mas and a median of 
1.2 mas. 

3.4.5. Transformation errors 

It is important to notice that any error in deriving 
the motions for the reference stars only translates into a 
global uncertainty of the coordinate system (which could 
show up as an offset of the center of mass from 0/0 or 
a net motion of the coordinate system). It will however 
not affect the accuracy of individual data points in this 
system. Only the selection of reference stars and trans- 
formation errors contribute to the errors of the individual 
data points. We estimated them by performing all coor- 
dinate transformations not only once but also with sub- 
sets of the available reference stars. The standard devia- 
tion of the sample of obtained astrometric positions was 
then included in the astrometric error estimate. The typ- 
ical uncertainty introduced by the transformations was 
quite small, namely 23 /ias for the NACO data. This 
is consistent with the fact that « 100 stars have been 
used of which each can be determined with an accuracy 
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of w 200 /zas. For the SHARP data we found a value 
of 100 n&s, again consistent with the characteristic single 
position error of « 1 mas. 

3.4.6. Differential effects in the field of view 

At the sub- mas level, there is a multitude of differential 
effects over the field of view that can influence astromet- 
ric positions. The most prominent ones are relativistic 
light deflection in the gravitational field of the sun, light 
aberration due to Earth's motion or refraction in the 
atmosphere. Since our analysis is based on relative as- 
trometry, the absolute magnitudes of the effects do not 
matter. Only the differential effects over the field of view 
can contribute to the positional uncertainties. 

• Over a field of view of 20" the differential effects 
of aberration c an be describ ed by a global change 
of image scale (jLindegren fc Bastian ,2006). Since 
we fit the image scale for each epoch separately, 
the differential aberration is absorbed into the lin- 
ear terms of the transformation and thus is not 
affecting the astrometry. The size of the effect for 
a small field of view with a diameter / amounts 
to f X v/c X cos^ where ^' is the angle between 
the observation direction and the apex point. For 
/ w 10" and V w 30km/s this yields « 1 mas at 
most. 

• The light defiection can be approximated by 
4 mas X cotvIf/2 where is t he angle between ob- 
serva tion direction and Sun (jLindegren fc BastianI 
I2006f ). The differential effect over 20" will not ex- 
ceed 100 /^as as long as ^ > 3.6°, which is guaran- 
teed for all our data. 

• From the usual refraction formula R — 44" tanz 
(for a standard pressure of 740mbar at Paranal) 
we find a differential effect of 4 — 8 mas over 20" 
or 2 — 4 mas over the field in which we selected the 
reference stars. The effect will be a change in one 
direction (towards zenith) of the image scale. Since 
the effect is at most quadratic over the field, it will 
be absorbed completely into the first and second 
order terms of the transformations. Note that it is 
crucial to allow also for skew terms, i.e. it is not 
sufficient to use a shift, rotation and scale factor 
only in the linear terms, but the off-diagonal terms 
in the transformation matrix are also required. 

3.4.7. Unrecognized confusion 

One important contribution to the position errors is 
the fact that stars can be confused and that sometimes 
the confusion is not recognized. This problem is more 
severe for the SHARP data than for the NACO data due 
to the lower resolution. Of course we excluded positions 
for which we know that they are confused. However, 
unrecognized confusion cannot be dealt with by princi- 
ple. We therefore simply accept that these events hap- 
pen. This means in turn, that we expect to find a re- 
duced > 1 when trying to describe the motions with 
smooth functions. In addition we note that for a suffi- 
ciently large amount of data points unrecognized confu- 
sion events should only lower the precision but not the 
accuracy since no global bias is expected. Still, if a con- 
fusion event happens during an unfortunate part of the 



orbit (for example at an end point or during pericen- 
ter passage) a bias in the results of an orbit fit can be 
introduced. 

3.4.8. Gravitational lensing 

Gravitational lensing might affect the measured posi- 
tions. A quantitative analysis shows that the effects are 
very small except in unusual, exceptional geometric con- 
figurations. For a star at a distance z Rq sufficiently 
far behind Sgr A* the angle of deflection as measured 
from Earth is 

z AGMuBH 
Ro 

where b is the impact parameter. For the GC, this evalu- 
ates to ~ 20 lias X z/b, indicative of a very small astro- 
metric effect unless z/b:^ 1 This rough estimate is con- 
sistent with the rigorous tre atment of the problem from 
iNusser fc B roadhurstI (|2004f ). who show that in order to 
achieve a displacement of 1 mas, a star at z w 1000 AU 
needs to have b ^ 2 mas « 16 AU. In our data set, none of 
the stars get close to the regime that gravitational lens- 
ing actually becomes important. Therefore, we neglected 
the effect. 

3.4.9. Comparison of error estimates with noise 

We were able to check how well our error estimates 
agree with the intrinsic noise of the data. For this pur- 
pose we fitted all measured positions of the reference 
stars with simple quadratic functions. After exclusion of 
3cr-outliers, we have calculated the reduced for each 
reference star. The mean reduced for the NACO data 
is 2.0±0.7, while for the SHARP data we obtained values 
between 0.5 and 2.0 with a mean of 1.0. 
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Fig. 7. — Final distribution of total astrometric errors for our 
data. Left column: SHARP data, right column: NACO data. 
Top row: reference stars, bottom row: S-stars. The curves show 
empirical fits to the histograms in order to determine the respective 
characteristic error as the peak of the distribution. 



Since our data set consists of two subsets (SHARP and 
NACO), each covering roughly the same amount of time, 
the relative weight of the two subsets matters. Given that 
we seem to underestimate the errors for NACO a bit, 
while the SHARP errors seem consistent with the noise 
in the data, we decided to apply a global rescaling factor 
of r = 1.42 to all NACO data points. This procedure 
adjusts the relative weight between the two subsets. Still 
we expect a reduced > 1 when performing orbit fits 
due to unrecognized confusion events. 
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In Figure [7] we show the final error distributions (af- 
ter rescaling all NACO errors with the global factor) for 
the S-stars and reference stars, both for the NACO and 
SHARP data. The characteristic error for a reference 
star in the NACO data is 360 /las, in the SHARP data it 
amounts to 760 /xas. For the S-stars, the histogram of the 
NACO errors has a peak also around 325 /xas and a tail 
towards larger errors, essentially telling us that for bright 
S-stars the astrometry is as good as one could hope for 
(since it is equally good as for the reference stars). The 
tail is due to the fact that many of the S-stars are faint 
(hence the statistical error is severe) and probably also 
unrecognized confusion events affect the statistical er- 
ror since confusion can alter the shapes of the images 
of faint stars. In the SHARP data, the typical S-stars 
error is 2 mas and the lower end of the distribution at 
« 1 mas is consistent with what could be expected from 
the reference stars. 

3.5. S2 in 2002 

Our data set covers the pericenter passages of several 
stars. Particularly important to our analysis is the one of 
the star S2. The star is one of the brightest in the sam- 
ple and we observed a full orbit (see Figure [0)1 . In 2002 
S2 passed its pericenter, thus changing quickly in veloc- 
ity throughout a period of a few months. These data 
are particularly useful for constraining the potential of 
the MBH. However, as we will now discuss, the photom- 
etry of the star near pericenter-passage is puzzling and 
may indicate that the positional information is affected 
by a possible confusion event with another star. Fig- 
ure [8] shows a K-band PSF-photom etrically determined 
light curve for the star ()Rank Il2007[) . It is clear that S2 
was brighter in 2002 than in the following years. There 
are several reasons why a star could change its apparent 
brightness. 




2003 2004 2005 2006 2007 
year 

Fig. 8. — The K-band magnitude of S2 as function of time in 
the NACO data, determined by means of PSF photomery (black 
data). For comparison the star S8 is shown (red data). 



1. In 2002, S2 was positionally nearly coincident with 
Sgr A* and thus confused with the NIR counter- 
part of the MBH. Typically, Sgr A* is fainter than 
tok = 17 and thus the extra- light from Sgr A* 
in quiescence is not sufficient to explain the ob- 
served increase in brightness of S2. However, 
Sgr A* is known to exhibit flares that can reach 
a brightness level that could account for the ob- 
served increase in brig htness (jCenzel et al.ll2003al : 



iTrippe et aLll2007D . In that case we would expect 
to see intra-night variability of S2 in the 2002 data. 
Assuming conservatively that we can determine the 
relative flux of S2 to Atok = 0.1 in each frame 
and given the brightness of S2 (mx ~ 14) we esti- 
mate that we would have noticed any variations in 
Sgr A* that exceed rriK ~ 16.5. Since we did not 
observe any intra-night variability we exclude that 
fiares from Sgr A* significantly contributed to the 
increased brightness of S2 in 2002. 

2. Intrinsic variability of S2 might explain the ob- 
served light curve. However, it is unlikely to be 
the correct explanation, since it would be a big co- 
incidence that the brightening happens during the 
pericenter passage. Also an eclipsing binary seems 
unlikely given the slow variation. 

3. The star could change its properties during 
the pericen t er p assage. While tidal heating 
(|Alexander I l2005f) cannot plausibly change the 
temperature of a star within a few months, the in- 
teraction of S2 with some ambient medium does not 
seem ruled out. Such an encounter would primar- 
ily change the surface temperature of the star and 
therefore would act nearly instantaneously. Effec- 
tively the light curve would then be a direct trace 
of the density of the surrounding gas encountered 
along the orbital path of S2. However, energeti- 
cally, this scenario seems unlikely: Given the max- 
imum velocity of S2 at pericenter {v » 8000 km/s), 
the ra dius of the star (r = Hi?©, iMartins et al.l 
(l200l 'l and assuming that the kinetic energy of 
the gas that hits the geometric cross section of 
the star is converted to radiation, one can estimate 
the number density n necessary to produce the ob- 
served brightness increase of Ato^ k, 0.5. We ob- 
tained n « lO^-'^cm"'^, which is unrealistically high, 
and so we do not favor this scenario. 

4. iLoeb I ()2004D proposed that the stellar winds of 
early-type stars passing their pericenters close to 
the MBH could alter the accretion flow onto 
Sgr A*. Such an event would produce a change 
in the brightness of Sgr A* on the timescale of 
months, compatible with Figure [51 However, 
IMartins et al.l (|2008( ) showed that the mass loss 
rate of S2 is too low for this mechanism to work. 

5. The extinction could be locally smaller than the 
average value. For instance, Sgr A* could remove 
dust in the interstellar medium in its vicinity. This 
hypothesis can be tested in the future by observing 
other S-stars passing close to Sgr A* during the 
pericenters of their orbits. 

6. The brightness of S2 could be affected by dust in 
the accretion flow onto the MBH. The dust would 
be heated by S2 and account for the excess bright- 
ness, a proposal that was used by iGenzel et al.l 
()2003bf ) to explain the MIR excess of S2/Sgr A*. 

7. The star could be confused with another star. If 
S2 had been located very close to another star in 
projection, the true nature of this encounter could 
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remain undiscovered, but the observed brightness 
of S2 would be increased. 

Of the three viable explanations (5 to 7), the first would 
not lead to astrometric biases, the others however would 
displace S2 artificially. Given the importance of the 2002 
data, we decided not to discard it completely but to es- 
timate the astrometric error assuming a confusion event, 
given the measured increase in brightness. 
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Fig. 9. — Simulation of a confusion event. The contour lines show 
by which amount a mx = 14 source is displaced if it is confused 
with a second source that has certain magnitude and that is located 
in a given distance. The units are pixels, the simulation assumed 
simple Gaussian point spread functions that are sampled as it is 
the case for the NACO detector in K-band. The area to the top 
right can be excluded since a relatively bright source a few pixels 
apart from the primary would produce an elongated shaped image 
(which is not observed for S2 in 2002). The line denotes the limit 
at which the major axis is 30% larger than the minor axis. The 
horizontal lines indicate the brightnesses that a secondary source 
would have needed to push the S2 brightness up by the observed 
amount for the observed magnitudes at the dates indicated. For 
each date a mean deflection can be read from this plot. That value 
is used as astrometric error for S2 at the given date. 

For this purpose we simulated confusion events. We as- 
sumed simple Gaussian point spread functions and sam- 
pled them as they are sampled by the 13mas/pix scale 
of the NACO camera in K-band. By polluting a pri- 
mary source with a fainter secondary source we gener- 
ated a confused stellar image. This was then fit by a 
two-dimensional Gaussian and the displacement from the 
position of the primary source was determined. We var- 
ied brightness ratio and distance between the two sources 
systematically, yielding a displacement map (Figure [5]) . 
This map allows the determination of the possible range 
of displacements if the brightness of the secondary source 
is known. The range can be constrained further, since a 
bright secondary source in a few pixels distance will lead 
to very eccentric images that would be easily detected in 
the data. We excluded all points that would lead to a 
stellar image of which the major axis is more than 30% 
larger than the minor axis. Thus, from the measured 
S2 fluxes, the known, unconfused brightness of S2 and 
the roundness of the S2 images, we were able to con- 
strain the astrometric bias due to confusion. For each 



date we looked up in figure [9] the possible range of as- 
trometric displacements given the observed brightness of 
S2, essentially determining the profile along a horizon- 
tal line in the plot. The mean of this distribution was 
then considered as an additional 2D error to be added 
to the respective astrometric errors for that date. The 
such obtained error terms ranged between 2.37 mas and 
3.76 mas 



0.04 F 




-0.02 [- 

0.04 0.03 0.02 0.01 -0.01 -0.02 

R.A. ('■) 

Fig. 10.— The 2002 data of S2. The grey symbols show the 
measured positions, the errors are as obtained from the standard 
analysis and are not yet enlarged by the procedure described in 
section 13.51 The black dots are the positions predicted for the 
observation dates using an orbit fit obtained from all data other 
than 2002. The blue shaded areas indicate the uncertainties in the 
predicted positions resulting from the uncertainties of the orbital 
elements and of the potential, taking into account parameter cor- 
relations. The little ellipse close to the origin denotes the position 
of the fitted mass and the uncertainty in it. This plot shows that 
the S2 positions are dragged for most of the data by Ri 10 mas to 
the NE; they are not biased towards Sgr A*. 



We checked whether the residuals of the 2002 data, rel- 
ative to an orbit fit to the data other than 2002, show 
some systematic trend (figure [10]) and found that in- 
deed all points appear to be shifted systematically by 
10 mas « 1 pix towards the NE. Still, this is hard to in- 
terpret. In particular, S2 does not appear systematically 
displaced towards Sgr A*. Extrapolating backwards the 
track of S19 that was observed from 2003 on shows that 
it also was located close to the S2 positions in 2002. 
Again, there is no indication that S2 would be displaced 
towards the extrapolated positions of S19. Also, S19 
with mK ~ 16.0 is too faint to account for the observed 
increase in brightness of S2. Any other star that poten- 
tially was close to S2 in 2002 (candidates are S23, S38, 
S40, S56) is even fainter. From figure [5] one can see that 
a star with ttik ~ 14.4—14.0 in a distance of 2 — 2.5 pixels 
would be required to account for the observed shift. Fur- 
thermore, that secondary source would have to move for 
a few months and for « 40 mas nearly parallel to S2. It 
is extremely unlikely that we have missed such an event. 

From this analysis, it is clear that the weight of the 
2002 data will influence the resulting orbit fits, since 
these points will systematically change the orbit figure 
at its pericenter. At the same time we have no plausible 
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explanation for the increase in brightness and the sys- 
tematic residuals in the 2002 data; in particular a confu- 
sion event seems unlikely. Thus, it is clear that using the 
2002 data will affect the results, but we cannot decide 
whether it biases towards the correct solution or away 
from it. Therefore we use in the following two options: 

a) we include the 2002 data with the increased error bars; 

b) we completely disregard the 2002 data of S2. 

4. ANALYSIS OF SPECTROSCOPIC DATA 

Most of the radial velocities were obtained with 
SINFONI. For the few non-SINF ONI data w e used 
the already p ublished values (jGhez et al.l l2003l : 
lEisenhauer et al.|[2003i) . 

From the SINFONI cubes we determined spectra by 
manually selecting on- and off-pixels for each S-star and 
calculating the mean of the on-pixels minus the mean of 
the off-pixels. The spectra were then used to determine 
the radial velocities of the respective stars at the given 
epoch. We only used spectra in which we were able to vi- 
sually identify the stellar absorption lines without doubt. 
The most prominent features are the Br-7 line for early- 
type stars and the CO band heads for late-type stars. 

Both line profiles are non-trivial, possibly biasing the 
result when using a simple Gaussian profile to fit the line. 
The bias can be avoided by crosscorrelating the spectra 
with a template and determining the maximum of the 
crosscorrelation. 

For the CO band heads we use d a template spec- 
trum from iKleinmann fc Halll ()1986l ). We used the well- 
established tool 'fxcor' which is part of NO AO-package in 
iraf. We identified the following stars as late-type stars: 
SIO, S17, S21, S24, S25, S27, S30, S32, S34, S35, S38, 
S45, S68, S70, S73, S76, S84, S85, S88, S89, Sill. 

Also for the early-type stars one might be worried that 
radial velocity measurements are biased due to a com- 
plex line profile. In particular, Br-7 might be affected by 
nearby He lines. We tested this for the bright star S2, 
by generating a template from our 2004 - 2006 data^: 
we estimated for all S2-spectra the velocities by simple 
Gaussian fits to the Br-7 lins. We then Doppler-shifted 
all spectra to the 0-velocity (using the iraf task 'dopcor') 
and coadded them (using the iraf task 'scombine'). This 
resulted in a first template for S2. With this template 
we crosscorrelated all individual S2-spectra in the wave- 
length range 2.08 — 2.20 /im (using the iraf task 'fxcor') 
and obtained better estimates for the velocities. With 
these new velocities we reassembled the template spec- 
trum. We stopped after this first iteration since the ve- 
locity differences had already converged to a mean de- 
viation of 0.2km/s with a standard deviation of 2km/s. 
This template spectrum is shown in Figure [TT] We used 
it to determine the final S2-velocities. Comparing the re- 
sults to the initial estimates of the velocities showed that 
the Gaussian fits were not notably biased. The mean ve- 
locity difference was 8km/s with a standard deviation of 
27km/s. Therefore we simply used the Gaussian fits to 
the Br-7 liufi for the other early-type stars. We identified 
the following stars as early-type: SI, S2, S4, S5, S6, S7, 
S8, S9, Sll, S12, S13, S14, S18, S19, S20, S22, S26, S31, 
S33, S37, S52, S54, S65, S66, S67, S71, S72, S83, S86, 

* The combined S2 s pectrum created in t his context was also 
the basis for the work of lMartins et al.l II2008I '). 
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Fig. 11.— The combined S2 spectrum from the 2004 - 2006 
SINFONI data, used as velocity template. 



Before the measured velocities can be used in a fit they 
have to get referred to a common reference frame. The 
most suitable choice is the LSR. We used standard tools 
to determine the corrections which for our data only de- 
pend on the observing date and the source location. The 
observatory's position on Earth does not matter at the 
level of 15km/s accuracy, since it leads to a correction 
< 0.5 km/s. 

4.1. Radial velocity errors 

All radial velocities crucially depend on an exact wave- 
length calibration. The errors in the radial velocities 
were estimated from the following terms: 

• The formal fit error. For radial velocities which 
were obtained from a cross correlation with a tem- 
plate spectrum, the formal fit error is given by the 
fit error of the peak in the cross correlation, which 
is calculated routinely with the cross correlation 
routine. For the data for which we fitted a simple 
line profile to the spectrum the formal fit error is 
also an output of the fit routine. The magnitude 
of this error depends on the spectral type and the 
SNR in the spectrum. For a bright late-type star, 
e.g. S35 with tuk ~ 13.3, the formal fit error can 
be as small as 10 km/s, for a bright early-type star, 
e.g. S2 with rriK ~ 14.0, a typical value is 30 km/s. 

• Accuracy of wavelength calibration for Br-7. 

We used the non-sky-subtracted data cubes in or- 
der to determine the positions of atmospheric OH- 
lines. Comparing those to the nominal positions 
allowed us to estimate the accuracy of the wave- 
length calibration in the range of Br-7 the He- 
lines around 2.11 /im. The rms of the OH-line posi- 
tions around their nominal positions yielded errors 
in the order of 2 — 3 km/s. 

• Accuracy of wavelength calibration for CO 
band heads. Since there are no OH emission lines 
at wavelengths longer than 2.25 /im, we used atmo- 
spheric absorption features in the non-atmosphcrc- 
divided spectra of the respective standard stars in 
order to asses the accuracy of the wavelength cali- 
bration at the wavelengths of the CO band heads. 
This was possible since our standard stars were 
early-type stars (spectral type around B5) that do 
not show spectral features at the region of interest. 
We divided the region from 2.25 /im to 2.40 /im into 
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short windows of AA = 0.05 /Ltm and cross corre- 
lated each with a respective theoretical spectrum 
of the atmosphere. The typical resulting deviation 
was measured to be lOkm/s. The accuracy of the 
procedure is limited however by the accuracy by 
which the individual deviations can be measured, 
which yielded a value of lOkm/s, too. So probably 
the calibration is even more accurate than lOkm/s 
and consistent with what is found for the accuracy 
of the calibration for the shorter wavelengths. 

• Uncertainty of the underlying spectrum. The 

GC region is highly confused. Therefore we did not 
use an automated procedure to extract the spec- 
tra from the data cubes but selected the respective 
signal and off pixels manually. Since there is no 
clear prescription for what the optimum way for 
that procedure is, we extracted each spectrum sev- 
eral times. This allowed us to estimate the error 
due to the selection of signal and off pixels. While 
for bright stars {mx ~ 14) this error term is below 
lOkm/s, it becomes dominant for fainter stars. For 
an early-type star of rriK ~ 15.5 a value of 100 km/s 
is common. 

Since the wavelength calibration is determined indepen- 
dently for all data sets, these errors will average out with 
an increasing database. 

5. ORBITAL FITTING 

The aim of the orbital fitting is to infer the orbits of 
the individual stars as well as information on the grav- 
itational potential. A Keplerian orbit can be described 
by the six parameters semi major axis a, eccentricity e, 
inclination i, angle of the line of nodes Q, angle from 
ascending node to pericenter uj and the time of the peri- 
center passage tp. If the orbit is only approximately 
Keplerian, these parameters should be interpreted as the 
osculating orbital parameters. The parameters describ- 
ing a simple point mass potential are the distance to the 
GC, Rq, the mass of the central object, Mmbh, its po- 
sition and velocity. Note that the potential might also 
be more complicated, for example due to an extended 
mass component or due to the corrections arising from 
the Schwarzschild metric. These parameters can be in- 
ferred from our data by orbital fitting. 

After 16 years of high-precision astrometry of the in- 
nermost stars in our galaxy and a few years of Doppler- 
based radial velocity measurements the accuracy of the 
available data has reached a level at which one might 
hope to detect deviations from the Keplerian orbits on 
which the stars apparently move due to the existence of 
the MBH at the dynamical center of the Milky Way. Such 
deviations may be due to relativistic effects or are the ef- 
fect of an extended mass component possibly residing 
in the vicinity of the MBH. Both cases are scientifically 
highly interesting. In order to analyze these effects we 
implemented a general orbital fitting routine that per- 
mits the fitting of orbits in an arbitrary potential and 
that can take into account also relativistic effects. 

For a 1/r potential it is well-known that the solutions 
of the equations of motion of test particles are (Kepler) 
ellipses. Assuming such a potential, orbits can be fit- 
ted by adjusting the orbital elements, since there is a 



straightforward prescription for the calculation of the po- 
sition and velocity vectors at any given time from the 
orbital elements. However, a more general approach is 
needed if an arbitrary potential determines the dynam- 
ics. Then the trajectory has to be determined numeri- 
cally. The problem can be described by the initial condi- 
tions of each test particle plus the parameters describing 
the potential. For each set of parameters a with re- 
spect to the measured data can be calculated. One seeks 
the parameter values which minimize the x^- This is 
a computationally demanding problem as at each step 
of the high-dimensional minimization the equations of 
motion are solved numerically. We chose t he high-level 
tool Mathematica (jWolfram ResearchI |2005[ ) for the im- 
plementation and tested it thoroughly, e.g. by comparing 
results with results obtained from the former routine that 
explicitly uses ellipse-shaped orb i ts an d that was used 
for the work of lEisenhauer et aTl (|2005f ). Some features 
of the new routine are: 

• The NIR fiares of SgrA* are believed to appear 
at the position of th e center of mass for the orbits 
(jGenzel et al.ll2003al ) . When a flare occurs it there- 
fore is reasonable to take the measured position of 
the flare into account and to identify it with the 
center of mass. This can be achieved by letting 
this measurement contribute to the of the flt. 
In total we measured 22 times a position of Sgr A* 
(at various brightness levels, typically at ttzk ~15). 
Note that with such a fit, while possibly constrain- 
ing the potential parameters better, one gives up 
the possibility of testing whether the center of mass 
and the NIR counterpart of Sgr A* coincide. 

• We implemented four relativistic effects: 

a) the geometric retardation due to the finite speed 
of light, also called the Roemer effect. This in- 
volves numerically solving the retardation equation 
^obs = — z(tcm)/c, whcrc z is the coordinate 
along the line of sight, in order to know the position 
and velocity of the star at the time of emission. 

b) the relativistic Doppler formula, giving rise to 
the so-called transverse Doppler effect, affecting 
only the radial velocities. 

c) the gravitational rcdshift due to the potential 
of the central point mass, altering t he conversion 
of line positions to radial velocities. iZucker et al.l 
(|2006D show that effects a) - c) might become vis- 
ible in the radial velocity measurements during a 
close periastron passage of a star. 

d) the first general relativistic correction to the 
Newtonian potential as given by the Schwarzschild 
metric: V{r) = -GMuBii/r + GMmbhI'^ /c^r^ 
where I is the orbital angular momentum of the 
star. 

Within the fitting routine all four effects can be 
turned on or off, or the strength of the effect can 
be used as a fit parameter where means the ef- 
fect is not present and 1 corresponds to the case in 
which the effect is as strong as expected from the 
theory. 

• We allow for additional mass components in the po- 
tential, described by an arbitrary number of addi- 
tional parameters, all of which can be either treated 
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as fixed or as free fit parameters. The additional 
mass components can be given either as a term in 
the potential or as a function describing the density 
as function of the spatial coordinates. In the latter 
case the routine determines the potential from the 
mass distribution by solving the Poisson equation 
V^V(r) = 47rGp(r). Here, one encounters either 
a case in which a closed solution for V{r) can be 
found or it might happen that for each set of pa- 
rameters for which is calculated during the fit 
the Poisson equation has to be solved numerically. 

• For some of the parameters of the problem there 
could exist independent measurements which one 
might want to take into account during the fit. An 
example is the position of the central mass. We 
used radio measurements of Sgr A* to determine 
the coordinate system and thus we expect the cen- 
tral point mass to reside in the origin of the chosen 
coordinate system. We therefore implemented the 
use of priors for any of the parameters, which can 
be done straightforwardly by including them into 
the calculation of x^- 

• Instead of fitting the semi-major axis, we fit the 
periastron distance p. This has the advantage that 
we can allow values of e < 0, effectively exchanging 
the role of major and minor semi axis. By using 
p the parameter space is compact and the fitting 
routine can smoothly pass e = 0. 

We followed the us ual approach when calculating the sta- 
tistical fit errors (|Pressl[l993) . For the given best fit 
solution at a certain set of values {pi} for the parame- 
ters we determine the Hessian matrix from the curvature 
of the x^-surface: d^y^ /dpidpj. The formal fit errors 
are the diagonal elements of the inverse of that matrix. 
Note that still these are only formal, statistical fit errors. 
Possible systematic errors come in addition to them. Pa- 
rameter correlations are taken into account by the matrix 
inversion. All orbital elements for a given star are corre- 
lated with each other and with the potential parameters. 
However, the other matrix elements describing correla- 
tions between orbital elements of different stars can be 
set to 0. This refiects the test particle approach in which 
one star can only influence the fit result for another star 
via its influence on the potential. We explicitly use the 
test particle approach also when calculating x^ for more 
than one star. It allows one to use several CPUs in paral- 
lel since the contributions to from the individual stars 
are independent. 

6. RESULTS 

In order to predict the motion of a star in a given gravi- 
tational potential one has to know six phase space coordi- 
nates, e.g. its position and velocity at a given time. Since 
the radial position is not measurable for any of the S-stars 
and only for a few the radial velocity is measured, one 
needs additional dynamical quantities. As such one can 
use accelerations, either in the proper motion or in the 
radial velocity. Also higher order derivatives (e.g. da/dt) 
of the astrometric data can be used as additional dynam- 
ical measurables. If more than six dynamical quantities 
are measures, the star can be used to retrieve information 
about the potential. 



This section is organized as follows: First, we check 
by polynomial fits (going up to third order), for which 
stars we can expect to find orbital solutions and which 
stars can contribute in the determination of the potential. 
Then we determine the potential, yielding also the orbits 
of the stars used in this step. Finally, we determine the 
orbits of the remaining stars in the given potential. 

6.1. Polynomial fits 

For stars for which a significant part of the orbit is sam- 
pled, the astrometric data cannot be described by poly- 
nomial fits anymore. Most prominently, in our data set 
this is S2 of which our astrometric measurements cover 
more than one complete revolution. For all other stars 
we report the polynomial fits to the astrometric data in 
the table in appendix [C] We also give there polynomial 
fits to the radial velocity data of those stars for which 
we were able to determine orbits. The order of the poly- 
nomials in all cases was chosen such that the highest 
order term still differed significantly (at the 5-cr level) 
from 0. Significances were calculated after rescaling the 
errors such that the reduced of the respective fit was 
1, which is a conservative approach. 

Astrometrically, we found significant da/dt (requiring 
at least a 5-(t level) values for the stars SI, S4, S12, S13, 
S14, S17 and S31. Significant astrometric accelerations 
(at the 5-CT level or above) were found in addition for 
S5, S6, S8, S9, S18, S19, S21, S23, S24, S27, S28, S29, 

533, S38, S39, S40, S48, S58, S66, S67, S71, S83, S87 
and Sill, where we checked that the acceleration vector 
actually points towards Sgr A*. 

We measured changes in the radial velocity for SI, S2, 
S4, S8, S13, S17, S19 and S24 (aU > 5-ct, except S24 with 
4.8cr). Additionally, we were able to determine radial 
velocities for S5, S6, S7, S9, SIO, SIX, S12, S14, S18, 
S20, S21, S22, S25, S26, S27, S29, S30, S31, S32, S33, 

534, S35, S37, S38, S45, S52, S54, S65, S66, S67, S70, 
S71, S72, S73, S76, S83, S84, S85, S86, S87, S88, S89, 
S92, S93, S95, S96, S97 and Sill. 

Summarizing, we expect 

• that the S2 data will dominate the problem of de- 
termining the gravitational potential; 

• that SI, S4, S8, S12, S13, S14, S17, S19, S24 and 
S31 can be used additionally to constrain the po- 
tential further; 

• that we can find orbits in addition for S5, S6, S9, 
S18, S21, S27, S29, S33, S38, S66, S67, S71, S83, 
S87 and Sill. 

The data for the stars for which we found orbital solu- 
tions is presented in figures [T2l and fT3l see also table [Cl 

6.2. Mass of and distance to Sgr A * 

Here and in the following we report always the fit re- 
sults including the (downweighted) 2002 data of S2 and 
excluding it. The coordinate system priors were used as 
given in equation 3] The fit errors reported are rescaled 
such that the reduced = 1. Note that these errors 
include the formal fit errors, taking into account param- 
eter correlations between the parameters reported here 
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Fig. 12. — The orbital data for the S-stars other than S2, the data of which is given in figure [T3l Left: The measured dechnations as 
function of time for the stars for which we were able to determine orbits together with the orbital solution. Middle: The same plot for right 
ascenscion. Right: The measured radial velocity for those stars for which we were able to measure changes in the radial velocity together 
with the orbital solutions. The radial velocities for the other stars are given in table ICl 



and the respective orbital elements determined simulta- 
neously. The systematic uncertainty due to the coordi- 
nate system is included here as well, since these param- 
eters were varied during the fits, too. The imp ortance of 
this was pointed out also by iNikiforovl (|2008[ ) . 

6.2.1. Ro and mass from S2 data only 

First, we used the S2 data only to determine a Kep- 
lerian gravitational potential (see figure [T3|) . Using the 
priors as obtained in equation 31 the fits yield the num- 
bers in the first and second row of table U) The two 
values for i?o differ by more than what the errors sug- 
gest; indicating that the 2002 data infiuences i?o- This 
confirms the presumption from section We exploited 
this further in figure [141 Assigning the 2002 data higher 
weights (smaller errors) pushes the distance estimate up, 
smaller weights lower it. 

Mass and distance are strongly correlated parameters, 
see Figure [151 The scaling of mass with Rq in our data 
set is a power law with A/mbh ^ Rq- For a purely 
astrometric data set one would have an exponent of 3 
and a complete degeneracy; the fact that the exponent 
is < 3 and that the degeneracy is not complete is due 
to the influence of the radial velocity information in our 
data set and due to the use of priors. The degeneracy 
can be understood qualitatively. Changing Rq effectively 
changes the conversion from measured angles (in mas) 
to physical lengths (in pc), i.e. changing i?o changes 
the semi major axis. Since the orbital period is well 
determined in our data, the mass has to change in order 
to fulfill Kepler's third law. 

The strong dependency means that the uncertainties 



for mass and distance are coupled. Fixing the dis- 
tance yields a very small fractional error on the mass 
of AA/mbh ~ 0.02AfMBH- This shows that the error of 
the fitted mass is completely dominated by the uncer- 
tainty in the distance. Once the distance is known, the 
mass immediately follows from the scaling relation 

Mmbh = (3.99±0.07| 

stat 

/ R \^ °^ 
X 10^ Mq —2- (incl.2002) , 
\ 8 kpc / 

Mmbh = (4.08 ± 0.09|,tat ± 0.39|rJ 

/ n \ 1-62 

X 10^ Mq — ^ (excl.2002) , (8) 

\ o iipC / 

where the error due to Rq corresponds to the fit error 
reported in table HI 

6.2.2. Position of the central point mass 

By construction the position of the radio source Sgr A* 
in our coordinate system is located at the origin. Since it 
is clear that Sgr A* is the MBH candidate we used this 
fact when applying the priors of equation |4l However, 
our data actually allows us to test this hypothesis. By 
leaving the position and proper motion of the mass com- 
pletely free, we can check how well the position of the 
mass coincides with Sgr A*. Using the S2 data only, no 
2D priors but the prior in Vz from equationjUwe obtained 
the numbers presented in the third and fourth row of ta- 
ble m We note that the mass is located within « 2 mas 
at the expected position. The current accuracy by which 
this statement hol ds is an improvement of a factor « 2 
over the work from lSchodel et all (|2003f) . 
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Fig. 13. — Top: The S2 orbital data plotted in the combined 
coordinate system and fitted with a Keplerian model in which the 
velocity of the central point mass and its position were free fit pa- 
rameters. The non-zero velocity of the central point mass is the 
reason why the orbit figure does not close exactly in the overlap 
region 1992/2008 close to apocenter. The fitted position of the 
central point mass is indicated by the elongated dot inside the or- 
bit near the origin; its shape is determined from the uncertainty 
in the position and the fitted velocity, which leads to the elonga- 
tion. Bottom: The measured radial velocities of S2 and the radial 
velocity as calculated from the orbit fit. 



We also report the S2-only fits when not using any 
coordinate system priors at all (rows 5 and 6 in table |4]) . 
This enlarges the errors on Rq and A/mbh substantially, 
the fit values however are not significantly different from 
the respective fits in which the u^-prior was applied. Not 
applying the u^-prior also shows a large uncertainty on 
Vz of « 50km/s; this parameter also is degenerate with 

Rq. 
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Fig. 14. — Fitted value of Rg for various scaling factors of the 
S2 2002 data, using a fit with the coordinate system priors. The 
factor by which the 2002 astrometric errors of the S2 data is scaled 
up strongly influences the distance. The mean factor determined 
in Figure [9] is fa 7, corresponding to Rq f» 8.1 kpc. 
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Fig. 15. — Contour plot of as function of Rq and central point 
mass. The two parameters are strongly correlated. The contours 
are generated from the S2 data including the 2002 data; fitting 
at each point all other parameters both of the potential and the 
orbital elements. The black dots indicate the position and errors of 
the best fit values of the mass for the respective distance; the blue 
line is a power law fit to these points; the corresponding function is 
given in the upper row of the text box. The central point is chosen 
at the best fitting distance. The red points and the red dashed 
line are the respective data and fit for the S2 data excluding the 
2002 data; the fit is reported in the lower row of the text box. 
The contour levels are drawn at confidence levels corresponding to 
Icr, 3(T, 5o", Tcr, 9cr. 



From the numbers it seems that the fit excluding the 
2002 data agrees better with the expectations for the co- 
ordinate system (equation 2]) than the fit including it. 
The latter is marginally consistent with the priors, while 
the former is fully consistent. This means that the 2002 
data not only affects Rq (which we want to measure and 
thus cannot use to judge the result) but also the position 
and velocity of the mass for which we have an indepen- 
dent measurement via the coordinate system definition. 
This argument points towards rejecting the 2002 data. 
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6.2.3. Position of the IR counterpart of Sgr A * 

At 22 epochs we have identified a source in the NACO 
data between 2003 and 2008 that might be associated 
with Sgr A*. In some cases, e.g. when a bright flare oc- 
curred, the identification seems unproblematic. In other 
cases, one cannot be sure that the emission is not due 
to an unrecognized star at or very close to the position 
of Sgr A*; an example is Figure [T] Due to this proba- 
bly very frequent confusion we expect that the measured 
positions are very noisy and we decided not to include 
them into the orbital fits. However, we checked whether 
the measured positions are compatible with the orbital 
fits. Fitting a linear motion model to the Sgr A* data 
we obtained 

a[mas] = (1.2 ± 0.8) + (0.15 ± 0.46) x (i[yr] - 2005.91) 
(5[mas] = (2.7 ± 0.7) - (0.73 ± 0.39) x (i[yr] - 2005. 9|^) 

The errors here are rescaled for a reduced of 1. The 
velocity errors are approximately a factor 5 larger than 
the priors from equation [H justifying our choice not to 
incorporate this data into the orbital fits. Given the un- 
certainties, the position of the IR counterpart of Sgr A* is 
consistent with the position of the central point mass. In- 
terestingly, that data seems to prefer a position of Sgr A* 
marginally North of the expected position, which is also 
the case for the orbit fits which include the 2002 data of 
S2. This weakens again the conclusion from section [6.2.2l 
that the 2002 data should be rejected. 

6.2.4. Ro and mass from a combined orbit fit 

Given the large uncertainties due to the 2002 data of 
S2, we decided to obtain more information about the 
potential by using a combined orbit fit and the coordi- 
nate system priors. For comparison, we also excluded 
S2 completely. We used the stars SI, S2, S8, S12, S13, 
S14. We selected these stars from the sample that can 
contribute to the potential (section [HT]) since for them a 
large fraction of the respective orbit is covered. We did 
not select S4 and S17 as they suffered confusion in the 
SHARP data. S19 was omitted because its time base is 
quite short still (the star was not detected before 2003). 
Since S24 would only contribute marginally to the po- 
tential, it was left out, too. Finally, we did not select 
S31, since the nearby sources S59 and S60 were confused 
with S31 in the earlier NACO data. N ot surprising, the 
final s ample contains the same stars as [Eisenhauer et aLl 
()2005f ) had reported orbits for. 

In order to balance the relative weights of the stars 
used, we had fitted the five additional stars first alone, 
leaving also the potential free (but applying the priors). 
While the such obtained flts were not of interest per se, 
they still provided a smooth, unbiased model for each 
star. Hence, we used the resulting reduced values 
to rescale the astrometric and radial velocity errors such 
that all stars yielded a value of 1. The scaling factors 
applied ranged from 1.20 to 2.33, the latter value being 
extreme and occurring for S13, which perhaps suffered 
from confusion in the SHARP data and of which the 
data in 2006/2007 was affected by confusion with S2. 
Our procedure guaranteed that such a star with a high 
astrometric noise would not contribute overly much to 
the combined x^- We obtained the results given in rows 
7, 8 and 9 of table |4| These numbers agree with each 



other within the uncertainties. The combined fit includ- 
ing the S2 2002 data also agrees with the corresponding 
S2-only fit. This is not true for the combined fit exclud- 
ing the S2 2002 data, which is hardly compatible with 
the respective S2-only fit. A possible reason is that the 
S2 data before 2002 is only relying on the SHARP mea- 
surements, which not only have larger formal errors but 
also is more affected by unrecognized confusion events 
than the NACO data. 

By fitting the combined data at various, fixed values 
of i?o we obtain again the scaling of mass and distance: 



A/mbh = (3.95 ± 0.06|stat ± 0.18|r„) 



xlO^Mn 



Ro 



2.19 



,8kpc 
Mmbh = (4.01±0.07| 

stat ± 0.18|rJ 



(incl.2002) , 



X 10*^ ( FTTI ) (excl.2002) , 



Skpc / 

Mmbh = (3.88 ± O.lOlstat ± 0.41|rJ 



X 10^ M(T 



8kpc 



3.07 



(excl.S2) , (10) 



6.2.5. Other systematic errors for Ro 

Beyond what was considered before, the physical 
model for the potential is another source of uncertainty. 
For example using a relativistic model instead of a Ke- 
plerian orbit model increased the distance by ARq = 
0.18 kpc (0.09 kpc) when including (excluding) the 2002 
data. This is consistent with the formal error on Rq. 
Since we do not detect explicitly relativistic effects, we 
stay with Keplerian orbits and consider the shift of the 
value as an uncertainty for Rq. Fitting a Plummer 
model (as in Section [^75]) instead of a point mass poten- 
tial increases the distance by a similar value: 0.14 kpc 
(0.03 kpc) when including (excluding) the 2002 data. 
The additional degree of freedom in this fit increased the 
formal uncertainty by 0.11 kpc added in squares. Finally, 
we adopted for the uncertainties of the potential an error 
of Ai?o = 0.25 kpc. 

An additional, systematic error is whether the use of 
priors (equation|4|) is correct. In order to address this, we 
repeated the combined orbit fits without the 2D priors. 
We obtained the numbers in rows 10 and 11 of table HI 
The influence of the priors on the value of Rq is relatively 
small (compare rows 7 and 8 with 10 and 11 in table |4]). 
We adopt for this source of uncertainty ARq = 0.10 kpc. 

Furthermore, rows 7, 8 and 9 of table [3] show that the 
uncertainty of the weights of the 2002 data from S2 in a 
combined flt alters Rq by ARq = 0.13 kpc. Deselecting 
S2 from the fits changes the result by ARq — 0.07 kpc. 
Finally, we assign ARo = 0.15 kpc for the uncertainties 
related to the selection of data. 

Adding up the uncertainties yields that the uncertainty 
of the distance to GC is still rather large with Atotai-Ro — 
0.35 kpc. Table [5] summarizes the error terms for Rq. 

6.2.6. Final estimate for Ro and mass 

We finally adopt the potential from the combined flt 
including the S2 2002 data, the difference to the one ex- 
cluding that data is negligible given the formal fit errors 
fsection [6.2.4p . This potential will be used in section [6^ 
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TABLE 4 

Results for the central potential from orbital fitting, from either S2 data only (rows 1-6) or a combined fit using in 
ADDITION SI, S8, S12, S13, S14 (rows 7 - 12). In rows 9 and 12, the combined fit was done without S2. The third column 

indicates whether the 2002 DATA FROM S2 WAS USED OR NOT; THE FOURTH COLUMN INFORMS ABOUT WHICH OF THE PRIORS FROM 

EQUATION [1] HAVE BEEN USED. 





Fit 


S2 


priors 


Ro 


.^'^MBH 


a 


S 


Va 


vs 


Vz 






2002 




(kpc) 


(1O«M0) 


(mas) 


(mas) 


(/jas/yr) 


(/ias/yr) 


(km/s) 


1 


S2 only 


yes 


2D, 


8.31 ± 0.33 


4.29 ± 0.35 


0.51 ± 0.64 


2.18 ± 0.89 


-5 ± 87 


119 ± 78 


0.8 ± 6.2 


2 


S2 only 


no 


2D, 


7.36 ± 0.43 


3.54 ± 0.35 


0.81 ± 0.66 


-0.63 ± 1.39 


-69 ± 91 


103 ± 81 


-0.8 ± 6.2 


3 


S2 only 


yes 


Vz 


8.48 ± 0.38 


4.45 ± 0.41 


0.37 ± 0.73 


2.33 ± 0.94 


76 ± 131 


231 ± 107 


0.8 ± 6.1 


4 


S2 only 


no 


Vz 


7.31 ± 0.45 


3.51 ± 0.36 


0.92 ± 0.75 


-0.84 ± 1.43 


-83 ± 137 


154 ± 114 


-0.9 ± 6.3 


5 


S2 only 


yes 


none 


8.80 ± 0.53 


4.93 ± 0.75 


0.31 ± 0.71 


2.44 ± 0.89 


74±127 


220 ± 107 


29 ± 36 


6 


S2 only 


no 


none 


6.63 ± 0.91 


2.85 ± 0.74 


0.96 ± 0.75 


-2.00±2.38 


-111 ± 148 


162 ± 115 


-42 ± 44 


7 


comb. 


yes 


2D, Vz 


8.33 ± 0.17 


4.31 ± 0.22 


0.80 ± 0.63 


2.19 ± 0.60 


-28 ± 71 


100 ± 68 


0.0 ± 5.0 


8 


comb. 


no 


2D, Vz 


8.20 ± 0.18 


4.22 ± 0.22 


1.07 ± 0.58 


1.54 ± 0.64 


-32 ± 73 


86 ± 71 


0.0 ± 5.1 


9 


w/o S2 




2D, Vz 


8.40 ± 0.29 


4.51 ± 0.49 


1.49 ± 0.99 


2.61 ± 1.37 


-66 ± 94 


-116 ± 94 


-1.3 ± 5.1 


10 


comb. 


yes 


Vz 


8.38 ± 0.16 


4.36 ± 0.21 


0.73 ± 0.65 


2.10 ± 0.61 


51 ± 106 


211 ± 97 


-0.4 ± 5.1 


11 


comb. 


no 


Vz 


8.22 ± 0.20 


4.25 ± 0.26 


1.22 ± 0.81 


1.59 ± 0.83 


± 133 


164 ± 123 


-0.5 ± 6.3 


12 


w/o S2 




Vz 


8.42 ± 0.31 


4.61 ± 0.55 


6.2 ± 2.0 


6.0 ± 1.9 


-335 ± 294 


-15 ± 281 


-1.0 ± 5.0 



TABLE 5 

Systematic errors for the distance to the GC, Rq- 



Error source 


Aflo(kpc) 


Fit error including position and velocity 




uncertainty of coordinate system 


0.17 


Assumed potential 


0.25 


Using priors or not 


0.10 


Selection of data 


0.15 


Total 


0.35 



Ro = 8.33 ± 0.17|stat ± 0.31 |,y, kpc . 



to determine the orbits of the other stars for which we 
expect to find an orbital solution. Hence, we find 

± 0.31 |,vs kpc . (11) 

It should be noted that this value is consistent within 
the errors with values published earlier (jEisenhauer et all 
[2003, 2005). The improvement of our current work is the 
more rigorous treatment of the systematic errors. Also 
it is worth noting that adding more stars did not change 
the distance much over the equivalent S2-only fit. For 
the mass we adopt 

MmBH = (3.95 ± 0.06|,tat ± 0.18|r„, etat ± 0.31 |r„, sys) 

2.19 



V8kpc^ 

: (4.31 ±0.38) X IO^A/q fori?o 



.33 kpc. (12) 



6.3. Testing for an extended mass component 

While Newtonian physics seems to describe the S-star 
system reasonably well, one actually expects to detect 
deviations from purely Keplerian orbits with accurate 
enough astrometric and spectroscopic data. There are 
two main reasons for this: 

• The relativistic effects as described in Sec- 
tion [5] lead to deviat i ons (iRubilar fc EckartI 120011 : 
iWeinberg et all [20051 : iGillessen et al.ll2006f ). Note 
that for S2 the pericenter advances by 0.18° per 
orbital revolution, not far from the precision of the 
orbit orientation in Table [71 

• In addition to the MBH a substantial amount of 
mass might reside in form of a cluster of dark 



stellar remnants around the M BH ([Morris 



Miralda-Escudc & Gould 2000; Munoetal 



Mouawad et al...2005;.Hopm an & Al exande"r[|2006f l 



1993 



2005 



This will also lead to a non-Keplerian orbit, with 
the pericenter precessing in retrograde fashion. 

Given our current data base S2 is the only star for which 
one can hope to find a deviation from a Keplerian orbit. 
Fitting a relativistic orbit to the S2 data yields a similar 
(158.5 compared to 158.7 for the Keplerian fit, both 
with 114 degrees of freedom). Allowing for an extended 
mass component in addition does not change much, 
typically we found ~ 157.4 (depending on the details 
of the model) at the cost of one additional free parameter. 

The simplest model for an extended mass component 
is a constant mass density p{r) described by 



Pir) = Po ■ 
More realistic is a power law model 



Pir) = Pa 



ro 



(13) 



(14) 



The power l a w mod el is motivated by the findings of 
iGenzel et al.l (|2003U ) who show that the stellar num- 
ber counts display such a density profile, which is also 
expected on theoretical grounds (jBahcall fc Wolj [19771 : 
lYoung|[T980( ). The parameters po and a are a character- 
istic density at the given radius and the power law index. 
We assum ed for the following a = —1. 4, a = —1.75 and 
a = -2.1 (Hopman fc Alexander 200 6|)- 

Similar investiga tions (Rubil ar fc Eckar"3 120011 : 
iMouawad et al.|[2005f ) used a Plummer model: 



p{r) 



3/iAfMBH 
core 



1 + 



-5/2 



(15) 



The free parameters of the Plummer model are the 
core radius fcoro and the mass parameter p, which cor- 
responds to the ratio of total extended mass versus 
mass of the central point mass. This model allows a 
convenient analytical description of the null hypothe- 
sis - no stellar cusp - and roughly describ es the surface 
light density distribution around Sgr A* (jScoville et all 
[200l iSchodel et al.l[200l . We adopt a core radius of 
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TABLE 6 

Results from S2 fits including an extended mass 
COMPONENT. The parameter. Tj describes the ratio of 
extended mass to the central point mass. The extended 

mass is accounted for in a spherical shell from the 
pericenter distance of S2 to the apocenter distance. The 
table shows the results for various potentials 



Fit 


Ro (kpc) r] 




incl. 2002 data of S2 


p — const 

power law, a — —1.4 
power law, a — —1.75 
power law, a — —2.1 
Plummer 


8.46 ± 0.25 0.029 ±0.026 
8.49 ± 0.26 0.020 ±0.017 
8.49 ±0.26 0.018 ±0.015 
8.52 ±0.27 0.015 ±0.013 

8.47 ± 0.26 0.025 ±0.022 




cxcl. 2002 data of S2 


p — const 

power law, a — —1.4 
power law, a — —1.75 
power law, a — —2.1 
Plummer 


8.00 ±0.33 0.018 ±0.028 
8.03 ±0.34 0.013 ±0.016 
8.03 ±0.34 0.012 ±0.014 
8.05 ±0.35 0.012 ±0.014 

8.01 ±0.33 0.016 ±0.023 



fcovc — 15mpc, which m atches the observed Hght pro- 
file (iMouawaid et al.ll2005f) . 

We fitted the S2 data for all three mass models and 
included in all cases the relativistic effects. The coordi- 
nate system priors were applied (equation 3]) and an ad- 
ditional prior was set on the Rq — 8.40 ± 0.29 kpc from 
the combined fit that excluded S2 completely (row 9, ta- 
ble [4|) . Any such fit can only test for mass inside the S2 
orbit; therefore we express the results in terms of mass 
enclosed between S2's apocenter (r = 0.230" = 8.9 mpc) 
and pericenter {r — 0.015" — 0.58 mpc) relative to the 
mass of the MBH and call this parameter i]: 

/•apo /• 

rj Mmbh — / dr dmn{r,m) (16) 

./peri J 

The results are shown in Table[6]from which we obtain 

= 0.021 ± 0.019|,tat ± O.OOeUod (incl. 2002) 
T]S2 = 0.014 ± 0.019|,tat ± O.OOSUod (excl. 2002). (17) 

The statistical fit error includes the uncertainties due to 
the coordinate system definition. The result corresponds 
to a 1-CT upper Hmit of 77 < 0.040 (0.033) and a 99% upper 
Hmit of < 0.066 (0.058) including (excluding) the S2 
2002 data, where the upper limits are defined such that 
the cumulated probability density function reaches the 
specified si gnificance level at the respective value for ry 
(jFeldman fc Cousins 1998). The (small) uncertainty in 
77 due to the model uncertainty has been excluded for the 
calculation of the upper limit since Table [6] shows that it 
affects rather the amplitude of ij than its significance. 

So the basic result of this study, improving m easure- 
ment u ncertainties b y a fa c tor of six over Schod el et al.l 
(|2002D : iGhez et all (20051) ; lEisenhauer eTall (|20"05l) . is 
that a single point mass potential is (still) the best de- 
scription of the data. Any deviations are smaller than a 
few percent of the point mass, within the orbits of the 
central S-star cluster. 

6.4. Stars with orbits 

Assuming the potential from section [6.2.6l we were able 
to determine orbits for the stars listed in section 16.11 
During these fits, each star was considered separately 
and the potential was fixed. This yielded a total of 26 



measured orbits as expected from section 16.11 An illus- 
tration of the (inner) stellar orbits is shown in Figure [TBI 
the orbital elements for all 26 stars for which we found 
orbits are summarized in Table [71 For the calculation 
of the errors quoted, all measurement errors (astrometry 
and radial velocities) were rescaled such that the reduced 
— 1- Furthermore, the uncertainties of the potential 
were included. 

As a double-check, w e ran Markov-Chain Monte Carlo 
(MCMC) simulations (|Tegmark et all 1200 4*) in order to 
asses the probability density distribution of the orbital 
elements in the six dimensional parameter space. Such 
a chain efficiently samples high-dimensional parameter 
spaces. The algorithm is simple: 

1. Choose a reasonable starting point in the parame- 
ter space and calculate for that point. 

2. Draw a random jump in the parameter space with 
the typical jump distance simultaneously for each 
parameter being the respective 1-cr uncertainty di- 
vided by the square root of the number of parame- 
ters (hence the mean jump distance corresponds to 
a 1-(T jump). The uncertainties are obtained from 
the Hessian matrix at the given point in parameter 
space. 

3. Calculate Xn f^'' the new point. 

4. If Xn < accept the new point, else accept the 
new point with a probability of exp(— (x^ — x^)/2). 

5. Store the new point if it is accepted in the buffer 
of the chain, otherwise store the old point. 

6. Go back to step 2. 

After running this chain for a many iterations, the distri- 
bution of points in the buffer of the chain measures the 
probability density distribution, which thus can be esti- 
mated by the chain. The actual implementation needs 
some extra tricks, e.g. for quicker convergence the pa- 
rameters should be chosen orthogonal to each other. In- 
terestingly, for a sufficiently long chain the result does 
not depend upon the chosen jump distance; that value 
influences rather how fast the chain samples the param- 
eter space. 

For each star we used the MCMC algorithm. Assum- 
ing some reasonable potential (e.g. as determined from 
a preliminary fit to the S2 data) we varied all six or- 
bital elements and checked whether the region in this 
six-dimensional parameter space which is reached by the 
chain is compact and reasonably well described by Gaus- 
sian functions (see Figure (TTj) . The advantage of doing 
so is mainly that, unlike a minimization routine that can 
be trapped in a local minimum, the MCMC simulations 
yield a global picture of the probability density distribu- 
tion. 

For all 26 stars for which we were able to determine 
an orbit the probability density distribution was well- 
behaved, i.e. in all cases the MCMC sampled a compact 
region in parameter space, the size of which was con- 
sistent with the expectation from the fit errors of the 
parameters. Examples are shown in Figure [TTl We con- 
clude that the orbital solutions presented in Table [7] are 
reliable. 
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0.4 0.2 0. -0.2 -0.4 

R.A. (") 



Fig. 16. — The stellar orbits of the stars in the central arcsecond for which we were able to determine orbits. In this illustrative figure, 
the coordinate system was chosen such that Sgr A* is at rest. 



Among the stars with orbital solution, six stars are 
late type (S17, S21, S24, S27, S38 and Sill). It is worth 
noting that for the first time we determine here the or- 
bits of late- type stars in close orbits around Sgr A*. In 
particular S17, S21 and S38 have small semi major axes 
of a « 0.25". The late-type star Sill is marginally un- 
bound to the MBH, a result of its large radial velocity 
(— 740km/s) at r = 1.48" which brings its total velocity 
up to a value w la above the local escape velocity. 

Furthermore we determined (preliminary) orbits for 
S96 (IRS16C) and S97 (IRS16SW), showing marginal ac- 
celerations (2.1(7 and 3.9a respectively). These stars are 
of special interest, since they were pro posed to mem- 
ber o f a clockwise rotating disk of stars (jPaumard et al.l 
[20M) . Similarly, we could not detect an acceleration 
for S95 (IRS16 NW). This excludes the star from being 
a me mber of the counter-clockwise disk (jPaumard et al.l 
[2006^, since in that case it should show an acceleration 
of ~ 150 /ias/yr'^, while we can place a safe upper limit 



of a < 30 /ias/yr^. 



300.8 
300.6 
300.4 
£ 300.2 
5 300.0 
299.8 
299.6 




1°) e fperi (yr-2100) 

Fig. 17. — Examples from the Markov-Chain Monte Carlo sim- 
ulations. Each panel shows a 2D cut through the six dimensional 
phase space of the orbital elements for the respective star. Left: 
Example of two well constrained and nearly uncorrelated param- 
eters. Middle: Example for two correlated parameters, which are 
nonetheless well constrained. Right: Example of badly constrained 
parameters, showing a non-compact configuration in parameter 
space. 
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TABLE 7 

Orbital parameters of those S-stars, for which we were able to determine orbits. The parameters were determined in the 
potential as obtained in section [6.2.61 the errors quoted in this table are the formal fit errors after rescaling them such 
that the reduced = 1 and including the uncertainties from the potential. the last three columns give the spectral 

TYPE ('e' for EARLY-TYPE STARS, 'l' FOR LATE-TYPE STARS), THE K-BAND MAGNITUDE AND THE GLOBAL RESCALING FACTOR FOR THAT 
STAR. Sill FORMALLY HAS A NEGATIVE SEMI MAJOR AXIS, INDICATIVE FOR A HYPERBOLIC ORBIT WITH e > 1. We ALS O CITE THE ORBITAL 
SOLUTIONS FOR THE STARS S96 AND S97 WHICH SHOWED ONLY MARGINAL ACCELERATIONS, SEE SECTION 17.3.11 



Star 




an 






e 






nn 




tp [yr-2000] 


T[yr] 


Sp 


rriK 


r 


SI 


0, 


.508 ± 


.028 


0, 


.496 ± 0, 


,028 


120.82 ± 0.46 


341.61 ± 0.51 


115.3 ± 2.5 


0.95 ± 0.27 


132 ± 11 


c 


14, 


,7 


1.49 


S2 


0, 


.123 ± 


.001 


0, 


.880 ± 0, 


,003 


135.25 ± 0.47 


225.39 ± 0.84 


63.56 ± 0.84 


2.32 ± 0.01 


15.8 ± 0.11 


c 


14, 


,0 


1.22 


S4 


0, 


.298 ± 


.019 


0, 


.406 ± 0, 


,022 


77.83 ± 0.32 


258.11 ± 0.30 


316.4 ± 2.9 


-25.6 ± 1.0 


59.5 ± 2.6 


c 


14, 


,4 


1.99 


S5 


0, 


.250 ± 


.042 


0, 


.842 ± 0, 


,017 


143.7 ± 4.7 


109 ± 10 


236.3 ± 8.2 


-16.4 ± 2.5 


45.7 ± 6.9 


c 


15, 


,2 


1.93 


S6 





.436 ± 


.153 


0, 


.886 ± 0, 


,026 


86.44 ± 0.59 


83.46 ± 0.69 


129.5 ± 3.1 


63 ± 21 


105 ± 34 


c 


15, 


,4 


1.45 


S8 





.411 ± 


.004 


0, 


.824 ± 0, 


,014 


74.01 ± 0.73 


315.90 ± 0.50 


345.2 ± 1.1 


-16.2 ± 0.4 


96.1 ± 1.6 


c 


14, 


,5 


1.20 


S9 





.293 ± 


.052 


0, 


.825 ± 0, 


,020 


81.00 ± 0.70 


147.58 ± 0.44 


225.2 ± 2.3 


-12.2 ± 2.1 


58 ± 9.5 


c 


15, 


,1 


2.23 


S12 


0, 


.308 ± 


.008 


0, 


.900 ± 0, 


,003 


31.61 ± 0.76 


240.4 ± 4.6 


308.8 ± 3.8 


-4.37 ± 0.03 


62.5 ± 2.3 


c 


15, 


,5 


1.54 


S13 


0, 


.297 ±0, 


.012 


0, 


.490 ± 0, 


,023 


25.5 ± 1.6 


73.1 ± 4.1 


248.2 ± 5.4 


4.90 ± 0.09 


59.2 ± 3.8 


c 


15, 


,8 


2.33 


S14 


0, 


.256 ±0, 


.010 


0, 


.963 ± 0, 


,006 


99.4 ± 1.0 


227.74 ± 0.70 


339.0 ± 1.6 


0.07 ±0.06 


47.3 ± 2.9 


c 


15, 


,7 


1.99 


S17 


0, 


.311 ± 0, 


.004 


0, 


.364 ±0, 


,015 


96.44 ± 0.18 


188.06 ± 0.32 


31945 ± 3.2 


-8.0 ± 0.3 


63.2 ± 2.0 


1 


15, 


,3 


2.46 


S18 


0, 


.265 ± 0, 


.080 


0, 


.759 ± 0, 


,052 


116.0 ± 2.7 


215.2 ± 3.6 


151.7 ± 2.9 


-4.0 ± 0.9 


50 ± 16 


c 


16, 


,7 


2.34 


S19 


0, 


.798 ± 0, 


.064 


0, 


.844 ±0, 


,062 


73.58 ± 0.61 


342.9 ± 1.2 


153.3 ± 3.0 


5.1 ± 0.22 


260 ± 31 


c 


16, 


,0 


2.31 


S21 


0, 


.213 ± 0, 


.041 


0, 


,784 ± 0, 


,028 


54.8 ± 2.7 


252.7 ± 4.2 


182.6 ± 8.2 


28.1 ± 5.5 


35.8 ± 6.9 


1 


16, 


,9 


1.55 


S24 


1 


.060 ± 0, 


.178 


0, 


.933 ± 0, 


,010 


106.30 ± 0.93 


4.2 ± 1.3 


291.5 ± 1.5 


24.9 ± 5.5 


398 ± 73 


1 


15, 


,6 


1.78 


S27 


0, 


.454 ± 0, 


.078 


0, 


.952 ± 0, 


,006 


92.91 ± 0.73 


191.90 ±0.92 


308.2 ± 1.8 


59.7 ± 9.9 


112 ± 18 


1 


15, 


,6 


1.79 


S29 


0, 


.397 ± 0, 


.335 


0, 


.916 ± 0, 


,048 


122 ± 11 


157.2 ± 2.5 


343.3 ± 5.7 


21 ± 18 


91 ± 79 


c 


16, 


,7 


1.92 


S31 


0, 


.298 ± 0, 


.044 


0, 


,934 ± 0, 


,007 


153.8 ± 5.8 


103 ± 11 


314 ± 10 


13.8 ± 2.2 


59.4 ± 9.2 


c 


15, 


,7 


1.97 


S33 


0, 


.410 ± 0, 


.088 


0, 


,731 ± 0, 


,039 


42.9 ± 4.5 


82.9 ± 5.9 


328.1 ± 4.5 


-32.1 ± 6.5 


96 ± 21 


c 


16, 


,0 


2.02 


S38 


0, 


.139 ± 0, 


.041 


0, 


,802 ± 0, 


,041 


166 ± 22 


286 ± 68 


203 ± 68 


3.0 ± 0.2 


18.9 ± 5.8 


1 


17, 


,0 


2.13 


S66 


1 


.210 ± 0, 


.126 


0, 


,178 ± 0, 


,039 


135.4 ± 2.6 


96.8 ± 2.9 


106 ± 6.3 


-218 ± 23 


486 ± 41 


c 


14, 


,8 


1.15 


S67 


1 


.095 ± 0, 


.102 


0, 


,368 ± 0, 


,041 


139.9 ± 2.3 


106.0 ± 6.1 


215.2 ± 4.8 


-305 ± 16 


419 ± 19 


e 


12, 


,1 


1.53 


S71 


1 


.061 ± 0, 


.765 


0, 


,844 ±0, 


,075 


76.3 ± 3.6 


34.6 ± 1.5 


331.4 ± 7.1 


-354 ± 251 


399 ± 283 


c 


16, 


,1 


2.44 


S83 


2, 


.785 ± 0, 


.234 


0, 


,657 ±0, 


,096 


123.8 ± 1.3 


73.6 ± 2.1 


197.2 ± 3.5 


61 ± 25 


1700 ± 205 


c 


13, 


,6 


1.23 


S87 


1 


.260 ± 0, 


.161 


0, 


,423 ± 0, 


,036 


142.7 ± 4.4 


109.9 ± 2.9 


41.5 ± 3.7 


-353 ± 38 


516 ± 44 


c 


13, 


,6 


0.94 


Sill 




-10.5 ± 7.1 


1 


,105 ± 0, 


,094 


103.1 ± 2.0 


52.8 ± 5.4 


131 ± 14 


-55 ± 70 




1 


13, 


,8 


0.94 


S96 
S97 


1.545 ± 0.209 
2.186 ± 0.844 


0.131 ± 0.054 
0.302 ± 0.308 


126.8 ± 2.4 
114.6 ± 5.0 


115.78 ± 1.93 
107.72 ± 3.15 


231.0 ± 9.0 
38 ± 52 


-376 ± 34 
175 ± 88 


701 ± 81 
1180 ± 688 


c 
c 


10.0 
10.3 


1.40 
1.15 



7.1. The distance to the Galactic Center 

Our estimate Rq = 8.33 ± 0.17|stat ± 0.31|syskpc 
equation HID is comp atible with our earlier work 
Eisenhauer et al.l 120031 l2005f) . While the underlying 
data base is partially identical, this work mainly im- 
proved the understanding of the systematic uncertain- 
ties. In particular, the astrometric data during the peri- 
center passage of S2 is hard to understand. This is an 
unfortunate situation, since that data potentially is most 
constraining for the potential. During the passage the 
star sampled a wide range of distances from the MBH, 
corresponding to a radially dependent measurement of 
the gravitational force acting on it. Probably only fu- 
ture measurements of either S2 or other stars passing 
close to Sgr A* will allow one to answer the question, 
whether the confusion problem close to Sgr A* is generic 
or whether 2002 was a unlucky coincidence. 

Besides stellar orbits, there are other techniques to de- 
termine Rq. a cla ssical one is to us e the distribution 
of globular clusters. iBica et all (l2006f ) applied this tech- 
nique to a sample of 153 globular clusters and obtained 
Rq = 7.2 ± 0.3. This value is only marginally compatible 
with o ur result. However, the error quoted bv lBica et all 
(|2006f ) corresponds to the formal fit error derived from 
their figure 4. Therefore, one might suspect that system- 
atic problems owed to the method were not yet included 
in the error estimate. 

The fact that the absolute magnitudes of red clump 
stars is known and that the red clump can be identi- 



fied in the luminosity function obtained from the ap- 
parent magnitudes of sta r s in the galactic bulge was 
used bv iNishiuvama et al.l (|2006[ ). These authors obtain 
Rq = 7.52 ± 0.10|stat ± 0.35|syskpc, where the statistical 
error is owed mainly to the uncertainty of the local red 
clump stars luminosities and the systematic error terms 
includes uncertainties in the extinction and population 
corrections, the zero point of photometry, and the fitting 
of the luminosity function of the red clump stars. This 
result is in agreement with our measurement, given the 
errors of both results. 

The known absolute magnitudes from RR Lyrae 
stars and Cepheids are the k ey to the work from 
iGroenewegen. Udalski fc Bond (|2008f ). Their result 
i?o = 7.94 ± 0.37|stat ± 0.26|statkpc is fully consistent 
with our result. The statistical error here is due to 
the photometric measurement errors, the zero point of 
photometry and the uncertainty of extinction correction. 
The systemtatic error includes the calibration of period- 
luminosity relations used and the selection effect, which 
could affect the result since only 39 Cepheids and 37 RR 
Lyrae stars have been used for this statistical approach. 

7.2. Limits on the binarity of Sgr A * 

It is interesting to see how our data limits the possi- 
ble existence of a second, intermediate mass black hole 
(IMBH) in the GC. Here, we do not aim at a rigorous 
treatment of the problem (which would be beyond the 
scope of this paper) but limit ourselves to estimates that 
appear reasonable given our findings. 
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The first constraint comes from the fact that the center 
of mass does not move fast. If the central mass were in 
orbit with an IMBH, the orbital reflex motion of Sgr A* 
might show up in our data. The upper limit on the ve- 
locity which we obtain from row 7 in tabic |4] corresponds 
to a line in a phase space plot of IMBH mass versus 
IMBH-MBH distance (Figure (TS]), separating configura- 
tions at smaller masses from systems with higher masses. 
From our data, we would not have been able to detect 
such an orbital motion of the MBH if the orbital pe- 
riod P were too short, namely much shorter than the 
orbital period of S2. We estimate that configurations 
with P > 5yr would be discoverable. Taken together, 
this excludes an area towards higher masses and larger 
distances. This constraint assumes implicitly that the 
stellar cluster rests relative to the MBH since it was de- 
rived in the combined coordinate system. Using the ve- 
locity calibration of the maser system would have yielded 
a slightly weaker constraint. However, an even stronger 
constraint comes from the radio measurements of Sgr A* 
(jRcid & Brunthaler 2004). The limit on the motion of 
radio Sgr A* in galactic latitude (vb = — 0.4± 0.9km/s) 
can also be used. Since this velocity limit is much smaller 
than the upper limit on the MBH motion from the stel- 
lar orbits, it is more constraining. Also for this data 
it seems reasonable to assume that only systems with 
P > 5yr would have been discovered. Similar arguments 
co nstraining the binarity of Sgr A* have been put forward 
bv [Hansen fc Milosavljevid ()2003[ ). whose results we also 
show in Figure [TSl 

Secondly, two black holes in close orbits will loose en- 
ergy via gravitational waves and thus spiral in. Demand- 
ing a hfe time of at least 10^ yr for the IMBH-MBH 
system excludes configurations towards smaller distances 
and higher masses. Dynamical stability can also be de- 
mand ed for the S-star cluster as such. Mikkola & McrriJ3 
(|2008l ) have shown that an IMBH with a mass of 
10~^ Mmbh in a distance of 1 mpc would make the S- 
stars cluster unstable. It is reasonable to assume that 
this also holds for larger masses and radii at least as 
large as the S-star cluster extends 1"). 

Based on simulations, Gua landris &: MerrittI ()2007D 
concluded that an IMBH will reach a stalling radius 
that is proportional to the mass of the IMBH: Ostaii = 
3.5 /ias X Mimbu[Mq] (for our values of mass and dis- 
tance). Since one does not expect an IMBH to reside at 
a much smaller radius, this puts another constraint on 
the IMBH-MBH binary. 

Finally, also the S2 orbit allows us to exclude part of 
the phase space. Motivated by the findings of section [^31 
and equation 1171 we simply assume that no mass larger 
than 0.02AfMBH can be hidden inside the S2 orbit. Ac- 
tually, also somewhat smaller radii than the pericenter 
distance rp of S2 are excluded, since this would still per- 
turb the orbit figure notably. We estimate that down to 
0.5rp no IMBH more massive than 0.02 Mmbh can reside. 

7.3. Properties of the stellar orbits 

We obtained orbits f or 20 early-type stars. This rel- 
atively large number - [ Eisenhauer et al.l (|2005D had six 
orbits. iGhez et al] (|2005D seven - allows us to assess dis- 
tributions of orbital parameters and study the proper- 
ties of the stellar orbits thereby characterizing the S-star 
population. 




10^ 10^ lO'* 10*^ 10' 

Mass 2"'^BH (M ) 



Fig. 18. — Constraints on the binarity of Sgr A* as function 
of the mass of the secondary black hole and the distance between 
the black holes. The shaded areas are excluded due to various ar- 
guments. The diagonal lines assume an orbital motion of Sgr A* 
around the secondary and correspond to velocity limits obtained 
either from the S-stars or the motion of Sgr A* ( Reid & Brunthale^ 
120041) . We estimate that only periods longer than 5 yr would lead 
to an observable effect, thus excluding an area towards higher 
masses and large distances. Demanding that the lifetime of the 
binary b lack hole excee ds 10^ yr yields another constraint (from 
[Hansen fc Milosavlievic] I I2003I) ). These authors also made similar 
arguments for the motion of the Sgr A*, the resulting constraint 
is replicated in this plot (denoted as HMOS). The stability of the 
S-star cluster puts a further constraint (Mikkola & Mcrritt 2008), 
as does the stalling radius found by gualandris & Merritt ( 2007), 
denoted as GM07. Finally also the S2 orbit excludes some part of 
the diagram, since it apparently is Keplerian. 



7.3.1. Orientations of orbital planes 

Figure [19] illustrates the orie ntations of the orbita l 
planes for all stars from Table [71 iPaumard "eiTall ([2006h 
suggested that the six stars S66, S67, S83, S87, S96 and 
S97 (E17, E15 (Sl-3), E16 (SO-15), E21, E20 (IRS16C) 
and E23 (IRS16SW) in their notation) are members of 
the clockwise disk. Our findings explicitly confirm this. 
All six stars have an angular distance to the disk be- 
tween 9° and 21° with a mean and standard deviation 
of 15° ± 4°. This is somewhat (a factor of 2) more than 
the disk thickness of 14° ± 4° found by IPaumard et al.l 
(2006). However, statistically the difference is not very 
significant and only the inner edge of the disk is sampled 
here. All six disk stars have a semi major axis of a w 1" 
and a small eccentri city (e » 0.2 — .4) in agreement with 
the estimates from IPaumard efall (j2006). The orbital 
plane of S 5 is also consistent with the disk given its dis- 
tance of 18°. However, the lower brightness {niK = 15.2) 
of the star and the higher eccentricity (e > 0.8) of the or- 
bit make it unlikely that S5 is a true disk member. The 
next closest star to the disk beyond the six disk stars 
and S5 is S31 with an angular distance of 27°. We also 
note that the orbital solutions for S96 and S97 derived 
from marginal accelerations are consistent with the disk 
hypothesis. Therefore we are confident in these orbits, 
too. 

We used a Rayleigh test (Wilkie! ll983l ) to check whether 
the distribution of orbital angular momenta for the 22 
other stars for which we found orbits is compatible with 
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Fig. 19. — Orientation of the orbital planes of those S-stars for which we were able to determine orbits. The orientation of the orbits in 
space is described by the orbital angular momentum vector, corresponding to a position in this all sky plot, in which the vertical dimension 
corresponds to the inclination i of the orbit and the horizontal dimension to the longitude of the ascending node Q. A star in a face-on, 
clockwise orbit relative to the line of sight, for instance, would be located at the top of the graph, while a star with an edge-on seen 
orbit would be located on the equator of the plot. The error ellipses correspond to the statistical 1 cr fit errors only, thus the area covered 
by each is 39% of the probability density function. Stars with an ambiguous inclination have been plotted at their more likely posit ion. 
The stars S66, S67, S83, S87, S96 and S97 which were suspected to be part of the clockwise stellar disk by IPaumard et al.l II2006I) at 
(n = 99°, i = 127°) actually are found very close to the positi on of the disk. The lat ter is marked by the thick black dot and the dashed 
lines, indicating a disk thickness of 14° ±4°, the value found bv IPaumard et aLl II2006I) . The orbits of the other stars ajre oriented randomly. 



a random distribution. We found a probability of ran- 
domness oi p — 0.74; meaning that the non-disk stars do 
not show a preferre d orbit orientation. Using the projec- 
tion m ethod from ICuesta-Albertos. Cuevas fc FraimanI 
(|2007t l we obtained p = 1.0. The same statement also 
holds when testing for randomness of the subset of early- 
type stars. 

7.3.2. Distribution of semi major axes 

Figure [20] shows the cumulative probability distribu- 
tion function (pdf ) for the semi major axes of stars which 
have semi major axis smaller than 0.5", thus excluding 
the stars that are identified to be members of the clock- 
wise disk. The statistic is limited still (15 stars make 
up this sample), but nevertheless the distribution allows 
us to estimate the functional behavior of the pdf n{a). 
Due to the small number of data points we did not bin 
the data but used a log- likelihood fit for n(a). We found 
n{a) ~ a'^'^^" '^. This can be conve rted to a numbe r 
density profile as a function of radius (jAlexander II2005D . 
We obtai n n(r) ^ ^-i i±o.3^ consistent with the mass 
profile in iGenzel et al.l (l2003bD who found p(r) ^ r~^-^ 
and with the newer work in TSchodel et all (l200l who 
found p{r) ^ r~^'^ for the innermost region of the cusp. 

7.3.3. Distribution of eccentricities 

The distribution of eccentricities allows us to estimate 
the velocity distribution. Figure [2T] shows the cumu- 
lative pdf for the eccentricities of those young (early- 
type) stars which are not associated with the clockwise 
stellar disk. Using again a log-likelihood fit, we find 
n(e) ~ g2.6±o.9^ r^Yie profile still is barely consistent with 
n(e) ~ e, corres ponding to an isotropic, thermal veloc - 
ity distribution (jSchodel et all l2003l : lAlexander I l2005f ). 




semi major axis (") 



Fig. 20. — The cumulative pdf for the semi major axis of the 
early-type stars with a < 0.5". The two curves correspond to the 
two ways to plot a cumulative pdf, with values ranging either from 
to (N-l)/N or from 1/N to 1. The distribution can be represented 
by n(a) ~a0-9±0-3. 



This would be the expectation for a relaxed stellar sys- 
tem. However, given that the maximal lifespan for B 
stars (< 10^ yr) is much shorter than the local two body 
relaxation (TBR) time (~ 10^ yr, Alexander (2005)) one 
does not expect a thermal distribution. In this light, it is 
interesting to notice that the distribution appears to be 
a bit steeper (i.e. peaked towards higher eccentricities) 
than a thermal distribution. This might be a first hint 
towards the formation scenario for the S-stars. For exam- 
ple, it is exactly what one expects in the bina ry capture 
scenario (|Perets. Hopman fc Alexander! 120071 ) . in which 
the S-stars are initially captured on very eccentric or- 
bits (e > 0.98), and then subsequent relaxation gradu- 
ally smears out the distribution of eccentricities towards 
a thermal distribution. From the time scales involved, 
one expects that the latter is not reached completely, 
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so a high eccentricity bias remains, which in turn might 
fit nicely together with our indication for a steeper than 
thermal eccentricity distribution. 

This lays out a very interesting perspective for the con- 
tinuation of the orbital monitoring. Increasing the statis- 
tics of the eccentricity distribution by determining more 
stellar orbits will allow us to test explicitly whether it 
truly deviates from a thermal distribution and thus pro- 
vides us with a quantitative test for formation scenarios 
of the S-stars. 




0.4 0.6 
eccentricity 

Fig. 21. — The cumulative pdf for the eccentricities of the early- 
type stars which are not identified as disk members. The two 
curves correspond to the two ways to plot a cumulative pdf, with 
values ranging either from to (N-l)/N or from 1/N to 1. The 
distribution is only marginally compatible with n{e) ^ e (dashed 
line), the best fit is n(e) ~ e2-6±0-9. 



7.4. Estimates of the extended mass component 

In addition to the population of stars not yet re- 
solved by current instrumentation a cluster of dark ob- 
jects - e.g. stellar mass black holes (SBHs) a ,s pro - 
posed inlMorrisI (1X9931) :lMiralda-Escude fc Gouldl (|2000D : 

is 



lMunoetal.1 (|2005f) : iHopman fc AlexandeH (|2006[) 
plausibly present in the GC. As shown in section 
the orbital data allows to test for such extended mass 
components. Here we investigate several theoretical and 
observational constraints on the extended mass distribu- 
tion and relate these to rj. We mostly assume that the 
extended mass distribution is due to SBHs with a m ass of 
Mi, = rUi^MQ with = 10 (jTimmes et al.lll996D . since 
this component is likely to make up most of the mass of 
a potential dark cluster (jAlexander Il2007i ). 

7.4.1. Stellar number counts 



iGenzel et afl ()2003b[ ) and ISchodel et all (|2007l ) have 
inferred a stellar density profile for the GC from 
completeness-corrected stellar number counts. Assum- 
ing that the luminous objects trace the total mass, the 
number density profile that is determined reliably on the 
> 0.01 pc scale can be extrapolated to the S2 orbit. We 
obtain 

3.7 X 10-4 X (^) , (18) 



'n 



This extrapolation is quite uncertain, since mass segrega- 
tion predicts that the SBHs should have a much steeper 
slope than the less-massive luminous stars in the cen- 
tral 0.01 pc, where the SBHs dominate the t otal mass 
(|Hopman fc Alexandeil [20061: 1 Alexander 1 12007| ) . There- 
fore both the mass-to-number ratio and the slope of the 



density profile are expected to have a significant radial 
dependence. 

7.4.2. T/ie drain limit 

The drain limit is a conservative theoretical upper 
limit of the number of compact objects that can exist 
in steady-state around a MBH. It is given by the con- 
dition that the number of SBHs that can be packed in- 
side any given radius in steady state has to be smaller 
than the number of SB Hs scattered into the MB H over 
the age of the Galaxy (jAlexander fc Livid [2004D . This 
can be translated into a theoretical limit for rj. Close to 
TOv, = 10 and using A/mbh = 4 x lO^M© and t = 10 Gyr 
the relation can be approximated by 



^ < 0.0011 X (^) 



-0.7 



(19) 



The drain limit could be violated for a non steady-state 
situation. Indeed, the existence of the young star disks 
with a relatively well-defined age of 6 My suggests that 
star formation in the GC is episodic. However, the 
amount of mass from SBHs would hardly exceed 10^ Mq 
even assuming an optimistic, top-heavy initial mass func- 
tion, given that the total amount of mass in the disks is 

7.4.3. Dynamical modeling of ttie dark cluster in the GC 

The expected degree of central concentration of SBHs 
around the MBH can be estimated by modeling the dy- 
namical evolution of a system with a, present-day mas s 
function similar to that of the GC (jAlexander 1 120051. 
Monte Carlo simulations of the GC using the Henon 
method and including also stellar collisions and tidal dis- 
ruptio ns (M. Freitag, priv. comm.; see also lFreitag et al.l 
(j2006[ )) but neglecting star formation yield a rather 
flat mass density profile of W{MqIy>c^) (r/0.01 pc)-^-^, 
which translates to ~ lO^"'. Due to the statistical 
nature of this method the density profile at the very cen- 
ter is not well determined. An alternative analytic solu- 
tion for the steady state distribution using a much more 
idealized fo rmulation of the mass segregation problem 
(jHopman fc Alexander ,2006) yields a similar result of 
77 ~ 5 X 10~*. However, in this method the fixed bound- 
ary conditions far from the MBH may artificially main- 
tain a high density in the center by preventing the ex- 
pansion of the system. Nevertheless, the fact that these 
two different methods yield similar results also consistent 
with the drain limit lends some credence to this estimate. 

7.4.4. Diffuse X-ray emission of a dark cluster 

A cluster of compact objects will accrete the surround- 
ing gas and thus lead to X-ray emission, which for current 
X-ray satellites (« 1") would be barely resolved. Indeed, 
the X-ray source at the position of Sgr A* is slightly ex- 
tended (jBag anoff et al. 2003) . We fit the radial profile 
of Sgr A* as reported by iBaganoff et al.l (j2003l ) by the 
superposition of a point source with a Gaussian width 
of CTpt = 0.375" (jBaganoff et al.ll2003f) and an extended 
component with a free width CToxt- We obtain as em- 
pirical descript i on fo r the surface brightness profile of 
IBaganoff et all (j2003f ). Figure 6: 

B{r) [cts/arcsec^] = 73.5 g-^''/^'"'' -H40.3 g-'''/^'^-' (20) 
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with CToxt = 1-05". Thus, we obtain for the extended 
luminosity (assuming the same spectral index of point 
like and extended component) ix,oxt = 1-95 x 10^'^erg/s, 
accounting for « 80% of the total X-ray luminosity. 

The expected X-ray luminosity of a single compact ob- 
ject is given by the mass accretion rate and the radiation 
efficiency. A simple esti mate is given by assuming Bondi 
accretion (iBondi 11195^ : 

Mb = 4TrX{GM^)^neHmpC-^ « lO^g/s , (21) 

where A = 1/4, rig = 26cm~^ the electron number 
density, fi = 0.7 the mean atomic weight, nip the pro- 
ton mass, Cs = \/5 fc Te/3 n the speed of sound and 
Te = 1.3 keV (jBaganoff et all [2001 . iPessah fc Melial 
(|2003f ) estimate the accretion rate by 

Mp^7:rl,pv^l0^g/s, (22) 

with the accretion radius race = 2GAf^,/w^g w 3 x 
10^^ cm. Using the density from above and the Kep- 
lerian velocity at r = 1" one obtains consistently Mp w 
109g/s«AfB. 

The radiation efficiency depends on the type of object 
considered (Hallcr et al. 1996). For neutron stars 10% is 
assumed (Pessah & Mclia 2003), since the accreted ma- 
terial will fall onto a hard surface and the energy re- 
leased can be radiated away, resulting in a luminosity of 
Lns ~ 10^^ erg/s. For SBHs due to the absence of a sur- 
face most of the emission will be ther mal bremsstrahlun g 
yielding only « 2 x lO^o erg/s (Haller et al."1996'). 
This shows that ix,cxt cannot be due to SBHs, since one 
would need 10^^ objects to explain the observed luminos- 
ity. In the case of neutron stars, one would need « 20000 
objects within r < 1" in order to account for the ob- 
served luminosity, corresponding to « 0.07. However, 
this num ber exceeds the estimate of th e segregated cusp 
model of iHopman &: Alexandeil (|2006D who predict only 
« 100 neutron stars there. 

7.4.5. X-ray transients in a dark cluster 

iMuno et al.l (|2005f ) report an overabundance of X-ray 
transients in the inner parsec of the GC compared to 
the overall distribution of X-ray sources. The sources 
are classified as X-ray binaries (XRBs). These authors 
suggest a dynamical origin of the XRBs, namely an ex- 
change of type Binary -I- SBH XRB + Star. The rate 
density for this reaction is 

7+ = rib S cTi (23) 

where nb is the density of binaries, 

E TTfl^ + 27raG(A/b + M^) /af (24) 

the exchange cross section and cti = (G' MMBH/3r)^/^ 
the ID velocity dispersion. According to iMuno et al.l 
()2005( ) the number of XRBs is limited by dynamical 
friction which yields a characteristic life time of r = 
10Gyr(M^/MQ)-i(r/pc)i/2. Refining this argument, 
we also take into account the back reaction XRB + Star — > 
Binary -I- SBH and assume for simplicity equal exchange 
cross sections for forward and backward reaction. Both 
effects together yield a rate density of 

= -nnxRB S(Ti H , (25) 

Z T 



where n is the number density of stars, rixRB the number 
density of XRBs.^ In equilibrium one has 7+ = 7_, 
which allows one to solve for nxRB- After integrating 
JT-XKB over volume out to 1 pc and assuming a = 0.1 AU, 
n = 10'5pc-3(r/pc)^^ rib = O.ln, Mb = 3Mq, NU = 
lOAf0 and = /^n one obtains the number of XRBs 
in the central parsec as iVxRB = 7 x 10^/*, where is 
the relative number of SBHs to ordinary stars. A certain 
fraction /x of those will shine up as X-ray transients: 
/xA^XRB = A'x- Calculating 7/ from this yields 

= — / AT:r^n{r)dr (26) 

A/mBH Jpori 

= 3.8xlO-MlM, (27) 

/X 

relating the number of X-ray transients in the central 
parsec with 77. Us ing the values fx ^ 0.01 and A^x = 4 

(|Muno et al1l2005| ) we obtain ri » 1.5 x lO"'^. 

A more detailed investigation bv lDeegan fc Navakshin I 
(2007) shows that within r < 0.7pc a cusp of SBHs with 
« 20000 members is consistent with the number of dis- 
crete X-ray sources in the GC. Converting this number 
for the assumed profile of n{r) ^ r^^/* into 77 yields 
rj « 2.1 X 10~^, which is very similar to our estimate in 
the previous paragraph. 

7.4.6. Further aspects 

There are at least three other aspects of an extended 
mass component in the GC which are worth mentioning 
but beyond the scope of this paper. 

• Star formation in the presence of a dark cluster. 
The process of star formation in the GC might 
be altered significantly by the presence of a sub- 
stantial dark component. The additional perturba- 
tive gravitational forces due to the SBHs might as- 
sist star formation since they increase the inhomo- 
geneities in a star forming gas cloud. On the other 
hand, close encounters between individual clumps 
and SBHs might result in disruption of the clumps. 

• Interaction of the spin of the MBH with the dark 
cluster. The spin of the MBH is subject to evo- 
lution by several processes. While gas accre- 
tion and major mergers can increase the spin, 
the accretion of SBHs tends to decrease th e spin 
(iHuehes fc Blandfor dI |200l; iGammie et all [2004) . 
assuming many random infall events of isotropi- 
cally distributed SBHs. Furthermore there is the 
general relativistic spin-orbit-coupling between a 
SBH and the MBH spin, leading to a change of 
the spin directi on of the MBH but not to a change 
of its modulus (jLodato fc Pringlell2006l ). 

• Dark matter. Dark matter, which is widely 
accepted in cosmology, might also show up in 
dynamic measurements in the GC. However, 
iGendin fc Primackl (|2004l ) show that the den- 
sity of the dark matter at 0.01 pc is pdm ~ 
3 X IQ^Mq/pc^, which is negligible compared 

^ The factor 1/2 takes care of the fact that for the back reaction 
either the SBH or the ordinary star of the XRB could be replaced 
by the interaction partner. 
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to the theoretically predic ted stellar density 

there (iHopman fc Alexandeil I2006D . See also 
iVasiliev fc Zelnikovl i200^ . 

7.4.7. Conclusions for an extended mass component 

The various estimates for 77 all consistently point to- 
wards an expected value of « 10""^ — 10""*, approximately 
two orders of magnitude smaller than what we can mea- 
sure with orbital dynamics today. Nevertheless, some 
astrophysical insights are possible. 

Among the most important scientific questions in the 
GC is the origin of the S-stars, being a population of ap- 
parently young stars close to the MBH ( Ghez et al.ll2003l : 
iMartins et al.l [20081 ) . One possible origin is that these 
stars have reached their current orbits by TBR. Then 
the S-stars would have an isotropic, thermal velocity dis- 
tribution, naturally explaining the observed random dis- 
tribution of angular momentum vectors (Figure fTOj) . The 
number of stars visible is by far too low to make TBR 
efficient enough to account for the present population 
of S-stars. A hypothetical cluster of SBHs could accel- 
erate the process. The C handrasekhar TBR timescale 
(|Binnev fc Tremaindll987b is given by 

^ 0-34 

'■^ G2(Af,)2n.lnA • ^^^^ 

For a power law cusp around a MBH, the velocity disper- 
sion and the density are related to each other. Assuming 
In A 10, a power law index of —3/2 (which is approxi- 
mately what is observed) and a population of stars with 
a single mass one obtains a relaxation time independent 
of radius 

« 1.8 X 105 yr 77-1 (^j , (29) 

Thus, if the S-stars formed at t he same epoch as the 
stellar disks 6 x 10^ years ago (jPaumard et al.l I2006D 
and reached their presen t-day orbits by TBR , one needs 
r] > 0.033 for = 10 (jTimmes et al.lll996[ ). This ex- 
ceeds the expectations by at least two orders of mag- 
nitudes. If the S-stars were not born in the presently 
observed disks, but in older, now-dispersed disks, one 
can use Equation with the typical age of B stars 
(« 5 X 10"^ yr). For = 10 this yields 77 > 3.5 x 10"^, 
which could be marginally compatible with the other es- 
timates for 77. 

In order to assess the expected progress we simu- 
lated future observations with existing instrumentation 
and similar sampling. Continuing the orbital monitor- 
ing for two more years will lower the statistical error to 
Aij « 0.01, corresponding to ~ 2 x 10^ yr. This means 
we will soon be able to test the hypothesis that the S- 
stars formed in the disks and reached their current orbits 
by TBR. Furthermore there is a chance to rule out any 
TBR origin of the S-stars observationally in the near fu- 
ture, namely when 77 < 3.5 x 10~^ is reached. 

8. SUMMARY 

We continued our long-term study of stellar orbits 
around the MBH in the Galactic Center. This work 
is based on our large, high quality data base which is 
based on high resolution imaging and spectroscopy from 
the years 1992 to 2008. The main results are 



• The best current coordinate reference system uses 
all available IR p ositions of the SiO maser stars 
(|Reid et all 120071) for the definition of the origin 
and assumes that the stellar cluster around Sgr A* 
is intrinsically at rest such that it can be used for 
the calibration of the coordinate system velocity. 
Having more measurements of the maser sources 
both in the radio and the IR domain we eventually 
will be able to directly tie the coordinate system 
velocity to radio Sgr A* with a sufficient precision. 
Then the intermediate step of cross calibration with 
the stellar cluster can be dropped and the coordi- 
nate system definition would be independent from 
the assumption that the stellar cluster is at rest 
with respect to Sgr A*. 

• We obtained orbits for 28 stars. Eleven of those 
can contribute to the determination of the gravi- 
tational potential, we used up to six. For the first 
time we were able to determine orbital parameters 
for six of the late-type stars in our sample. Fur- 
thermore ,_we_confirm unambiguously the earlier re- 
port ( Paumard et al.ll2b06.1 that six of the stars are 
members of the clockwise disk. 

• Overall, we improved measurement uncertainties 
by a factor of six over the most recent s et of Galac- 
tic Center papers ( Schodcl et al. 2002; iGhez etldl 
[2005; Eis enhauer et al.l 120051) . A single point mass 
potential continues to be the best fit to these im- 
proved data as well. The main contribution to the 
error in the mass of Sgr A* and the distance to 
the Galactic Center are systematic uncertainties. 
While the value of the mass is driven by the dis- 
tance estimate, the latter is subject to many sys- 
tematic uncertainties that amount to 0.31 kpc. The 
statistical error now decreased to 0.17 kpc and be- 
came smaller than the systematic one. The most 
fruitful way to overcome current limitations would 
probably be the observation of another close peri- 
center passage of an S-star. Our current best values 
are: 

M= (3.95 ± 0.06|,tat ± 0.18|r„, stat ± 0.31|r„, sys) 

= (4.31 ± 0.36) X 10^ Mq for i?o = 8.33 kpc 
Ro = 8.33 ± 0.17|,tat ± 0.31|,y, kpc (30) 

It should be noted that this value is consistent 
within the errors with valu es published earlier 
([Eisenhauer et all [20031 [2005D . The improvement 
of our current work is the more rigorous treatment 
of the systematic errors. Also it is worth noting 
that adding more stars did not change the distance 
much over the equivalent S2-only fit. 

• We have obtained an upper limit for the mass en- 
closed within the S2 orbit in units of the mass of 
the MBH: 

r/ = 0.021 ± 0.019|stat ± O.OOeUodoi • (31) 
which corresponds to a la upper limit of 77 < 0.040. 
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A. LIST OF NACO DATA SETS 



Date 


Band 


mas 
pix 


DIT 


NDIT 


4k 
TT 


4k 
TT 


Date 


Band 


mas 
pix 


DIT 


NDIT 


J* 

ir 


4t 
ir 








[sj 




frames 


S- stars 






isj 




frames 


S-stars 


onno OK 


i\.S 




u.o 


o 





35 




XI 


1 Q 






OD 


y / 


zUUz,o4 


Ks 


27 


20 


3 


20 


84 


zUUu.4i 


H 


13 


17 


2 


48 


80 


onno on 


Ks 


13 


15 


1 


1 1 


64 


o f 1 n in 


H 


13 


17 


2 


48 


53 


onno /I f 1 n 


Ks 


13 


15 


1 


86 


70 


zUUo.4;.t 


Ks 


13 


17 


2 


94 


58 


onno /I 1 o 


Ks 


13 


15 


1 


58 


47 


onnii? c ^ 

ZUU6.57 


H 


13 


60 


1 


32 


55 


onno ATA 


Ks 


13 


15 


1 


45 


83 


onniT c o 

ZUU6.58 


Ks 


13 


2.4 


14 


38 


75 


onno 


ts^s 


13 


15 


/I 

4: 


60 


83 


onnc c 


H 


13 


17 


2 


44 


57 


onno i^K 


TT 

rl 


13 


15 


/I 

4 


25 


86 


onnK TOK 
ZUUD. / ZD 


±vs 


13 


17.2 


2 


A s 


85 


onno 


Ks 


13 


15 


4 


20 


84 


onnc TOO 
zUUb. 1 zo 


Ks 


13 


17.2 


2 


40 


89 


onnQ oi 


H 


13 


20 


1 


32 


84 


onnc 
zUUb. ID 


Ks 


13 


17.2 


2 


24 


88 


onno on 


H 


13 


10 


6 


30 


89 


onnc TO 
zUUD. /8 


Ks 


13 


17.2 


2 


48 


93 


onno oc 


Ks 


13 


5 


12 


8 


82 


onnc on 
zUUb.oU 


Ks 


13 


2.4 


14 


48 


87 


onno A A 


Ks 


13 


5 


3 


102 


81 


onnT 1 T 

2007. 17 


Ks 


13 


12 


3 


32 


100 


zuuo.4oi 


TT 

rl 


13 


20 


1 


116 


90 




H 


13 


17 


2 


102 


105 




Ks 


1 


20 




182 


oD 


2007.21 


Ks 


13 


2.4 


14 


48 


100 


zUuo.404 


H 


13 


20 


1 


150 


85 


2007.214 


Ks 


13 


17.2 


2 


96 


100 


2003.454 


Ks 


13 


10 


2 


208 


81 


2007.252 


Ks 


13 


10 


3 


48 


99 


2003.55 


H 


13 


20 


3 


72 


81 


2007.255 


H 


13 


10 


3 


96 


103 


2003.676 


H 


13 


20 


3 


32 


81 


2007.255 


Ks 


13 


10 


3 


63 


97 


2003.678 


H 


13 


2 


30 


32 


65 


2007.46 


Ks 


13 


17.2 


2 


110 


100 


2003.76 


Ks 


13 


5 


12 


34 


89 


2007.54 


H 


13 


10 


3 


48 


96 


2004.24 


H 


13 


10 


3 


41 


85 


2007.55 


H 


13 


10 


3 


96 


106 


2004.33 


H 


13 


15 


2 


73 


92 


2007.69 


H 


13 


17 


2 


48 


104 


2004.35 


Ks 


13 


10 


3 


52 


70 


2007.69 


Ks 


13 


17.2 


2 


48 


100 


2004.44 


H 


13 


15 


2 


48 


94 


2007.692 


Ks 


13 


17.2 


2 


48 


100 


2004.51 


Ks 


13 


30 


1 


272 


86 


2008.15 


Ks 


13 


17.2 


2 


48 


101 


2004.52 


H 


13 


30 


1 


48 


82 


2008.20 


Ks 


13 


17.2 


2 


68 


106 


2004.57 


H 


13 


15 


2 


47 


70 


2008.27 


Ks 


13 


17.2 


2 


96 


93 


2004.57 


Ks 


13 


15 


2 


92 


46 


2008.46 


Ks 


13 


17.2 


2 


96 


101 


2004.66 


Ks 


13 


15 


2 


100 


89 


2008.47 


H 


13 


17.2 


2 


65 


88 


2004.73 


H 


13 


25 


1 


16 


92 


2008.60 


Ks 


13 


17.2 


2 


90 


104 


2005.27 


Ks 


13 


2 


15 


48 


95 
















2005.37 


Ks 


13 


2 


15 


71 


91 


DIT: single 


detector integration time 






2005.47 


Ks 


13 


10 


2 


77 


83 


NDIT: number of single integrations per image file 




2005.58 


Ks 


13 


15 


4 


23 


93 
















2005.67 


Ks 


27 


30 


1 


19 


54 
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B. POLYNOMIAL FITS TO THE S-STARS DATA 

The following table lists the polynomial fits to the S-stars data (except S2 which is not well described by polynomial 

fits) . For stars with a significant (at the 5-cr level) astromctric acceleration pointing towards Sgr A* wc report quadratic 
fits. For stars with significant da/dt we report the cubic fit. Otherwise linear fits are given. Similarly, for stars for 
which detected changes in the radial velocities, we report linear fits. For stars for which we determined orbits but did 
not detect changes in v^Tad we report the weighted averages. 



to [yr] for {a, S) 
to [yr] for Vz 


[masj — 
S [mas] — 
Vz [km/s] = 


SI, 14.7 
2000.41 
2005.77 


S7 2 + 4'i -1- n Q 7n + n 1 's^i A/ -I- (Sfi^ + n 012') ( /xt)'^ -\- ( —0 048'^ + n nnsfiV A/ 

(-125.4 ± 0.5) + (-32.04 ± 0.16)A< + (1.080 ± 0.023)(At)^ + (0.0398 ± 0.0032)(A0'' 
(1094.5 ± 8.4) + (-45.8 ± 4.2)At + (13.6 ± 2.6)(At)^ 


S4. 14.4 
2003.07 
2006.40 


(^269 9 + n -1- M 5 95 + 16~lAt -1- 544 + 070U A/^^ ^ ('—0 0260 + OlOlVA/^l'^ 
(124.8 ± 0.1) + (-0.36 ± 0.08)At + (-0.129 ± 0.023)(At)= + (0.0134 ± 0.0039)(At)3 
( — 687.9 ± 13.3) + ( — 66.7 ± 9.1)At 


S5, 15.2 
2005.05 
2006.40 


(352.0 ± 0.3) + (-4.93 ± 0.12)At + (-0.510 ± 0.078)(At)^ 
(200.4 ± 0.4) + (8.00 ± 0.15)At + (—0.341 ± 0.094)(Ai)^ 
129 ± 39 


S6, 15.4 
2005.05 
2006.40 


(484.6 ± 0.1) + (6.38 ± 0.04)At + (-0.128 ± 0.026)(Af)2 
(99.1 ± 0.2) + (0.74 ± 0.05)At + (—0.075 ± 0.031)(At)^ 

118 ± 21 


S7, 15.3 
2000.41 


(533.2 ± 0.2) + (-4.03 ± 0.04)At 
( — 29.8 ± 0.3) + (-3.08 ± 0.05)At 


S8, 14.5 
2000.41 
2005.77 


(336.1 ± 0.2) + (15.03 ± 0.08)At + (-0.347 ± 0.009)(At)=^ 
(-212.6 ± 0.2) + (-14.60 ± 0.10)At + (0.228 ± 0.010)(At)^ 
— (53.6 ± 7.6) + (—31.2 ± 5.3)At 


S9, 15.1 
2001.46 
2006.40 


(181.1 ± 0.3) + (1.65 ± 0.16)At + (-0.254 ± 0.021)(At)^ 
(-335.9 ± 0.3) + (-8.69 ± 0.15)At + (0.634 ± 0.019)(At)^ 
614 ± 27 


SIO, 14.1 
2001.46 


(64.2 ± 0.1) + (-5.14 ± 0.03)At 
(-381.9 ± 0.1) + (2.96 ± 0.03)At 


Sll 14.3 
2000.41 


(142.4 ± 0.3) + (8.79 ± 0.05)At 
(-548.9 ± 0.2) + (-5.39 ± 0.04)At 


S12, 15.5 

2002.07 

2005.77 


(-66.4 ± 0.3) + (4.73 ± 0.21)At + (1.066 ± 0.052)(At)^ + (-0.1002 ± 0.0066)(At)=* 
(252.0 ± 0.3) + (29.01 ± 0.15)At + (—1.606 ± 0.028)(At)^ + (0.0745 ± 0.0044)(At)^ 
318 ± 7 


S13 15.8 

2001.46 

2006.40 


( — 187.3 ± 1.6) + (15.22 ± 0.79)At + (3.954 ± 0.129)(At)^ + ( — 0.0880 ± 0.0206)(At)^ 
(-32.0 ± 1.0) + (44.79 ± 0.48)At + (-0.085 ± 0.068)(A«)2 + (-0.3469 ± 0.0116)(At)3 
(-7.3 ± 22.8) + (186.7 ± 14.7)At 


S14 15.7 

2000.41 

2006.40 


(64.2 it 2.4) + (17.48 ± 0.89)Ai + (2.838 ± 0.085)(At)'^ + ( — 0.2773 ± 0.0142) (At)'"* 
(36.7 ± 2.2) + (14.13 ± 0.81)At + (2.485 ± 0.076)(At)^ + (-0.2473 ± 0.0129)(At)^ 
300.3 ± 25.2 


S17, 15.3 

2001.46 

2005.85 


(10.3 ± 1.0) + (2.81 ± 0.47)At + (-0.373 ± 0.072)(At)^ + (0.1068 ± 0.0139)(Ai)'^ 
(-162.1 ± 1.0) + (20.04 ± 0.39)At + (1.408 ± 0.068)(At)^ + (—0.1411 ± 0.0146)(At)^ 
(594.7 ± 5.2) + (-84.0 ± 4.6)At 


S18, 16.7 

2005.47 

2006.40 


(-202.5 ± 0.5) + (-17.08 ± 0.18)At + (0.946 ± 0.115)(At)^ 
( — 70.2 ± 0.6) + ( — 18.22 ± 0.22)At + (0.008 ± 0.146)(At)^ 

-257 ± 53 


S19, 16.0 

2005.98 

2006.40 


(38.7 ± 1.1) + ( — 10.27 ± 0.31)At + ( — 2.461 ± 0.216)(At)^ 
(-117.3 ± 0.9) + (-19.78 ± 0.26)At + (6.068 ± 0.177)(At)^ 
(-2314.9 ± 36.4) + (10.4 ± 15.4)At + (88.0 ± 13.7)(At)^ 


S20, 15.7 
2005.98 


(220.8 ± 0.6) + (-4.94 ± 0.31)At 
(109.5 ± 0.4) + (-6.29 ± 0.21)At 


S21, 16.9 

2005.47 

2006.40 


(-334.1 ± 0.2) + (-11.19 ± 0.06)Af + (1.001 ± 0.039)(At)^ 
(-128.7 ± 0.4) + (4.15 ± 0.10)At + (0.598 ± 0.068)(At)^ 
410 ± 12 


S22, 16.6 
2005.47 


(191.2 ± 0.2) + (22.81 ± 0.11)At 
(-268.4 ± 0.3) + (-6.89 ± 0.14)At 


S23, 17.8 
2005.47 


(307.4 ± 1.1) + (-13.81 ± 0.34)At + (-0.953 ± 0.196)(At)^ 
(-89.1 ± 0.9) + (-11.17 ± 0.23)At + (0.525 ± 0.153)(At)^ 


S24, 15.6 

2001.46 
200().4U 


(-177.1 ± 0.3) + (6.32 ± 0.14)At + (0.065 ± 0.019)(At)^ 

(-506,7 + 0.3) J- fin, 40 + n.i3)A/ - ((^.2V^ + n.nmifA/i^ 

^-^S24.9 1 ().5) - (-23.0 ± 4.K)(A;,) 


S25, 15.2 
2001.46 


(-95.3 ± 0.3) + (-2.88 ± 0.06)Af 
(-426.4 ± 0.3) + (1.08 ± 0.05)At 


S26, 15.1 
2001.46 


(514.4 ± 0.2) + (6.39 ± 0.05)At 
(440.6 ± 0.2) + (0.81 ± 0.05)At 


S27, 15.6 

2001.46 

2005.77 


(146.9 ± 0.3) + (0.66 ± 0.15)At + (-0.032 ± 0.019)(At)^ 
(523.1 ± 0.3) + (3.74 ± 0.14)At + (-0.112 ± 0.018)(At)^ 
-114 ±3 
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Name, rriK 

f-i'rl frir f ,-\ A\ 

Cq [y 1 J i( )i { (\ . () ) 


a [mas] = 

[km/h] = 


S28 17 1 
2001.46 


(427.6 ± 0.6) + (12.09 ± 0.27)At + (-0.648 ± 0.036)(At)^ 


S29, 16.7 
2007.94 


(-206.5 ± 0.2) + (1.79 ± 0.16)At + (0.256 ± 0.036)(At)^ 

/'c;10'^-l-n'^^_L/' t a -l-n9s'iA/-i-^ n 7Q9 -I- n nfi/i \ i A 
\^oiy.o Hz yj.o) -t- ^ — ift.yo xj.j^o j ial ^ — u./y^i u.ud4: j j 

-273 ± 38 


2000.41 


(384.3 ± 0.1) + (3.37 ± 0.02)At 


S31, 15.7 
zuu J- . oy 
2006.40 


(-303.7 ± 0.5) + (5.89 ± 0.23)At + (0.561 ± 0.040)(Ai)^ + (0.0269 ± 0.0068)(At)^ 

/"ooi fi-|-0'^1-U(' fifi -4- 94^A+ -U r "^47 -1- 04'^^('A + ^■^ -4- ( 0491 -1- 00fi>4^('A^■^^ 

-366 ± 23 


OO.i, ID.D 

2005.43 


^ — OZO.O i U.lj -f- — 41:. if i U.UOJiAE 

(-356.2 ± 0.2) + (0.62 ± 0.07)At 


OOO, ID.U 

2000.41 
2006.42 


( /111 7-l-nc:^_L/■ ^ "i. -\- r\ or\\ A + A- 1 r\ 0^"^ A- r\ r\or\\ 1 

\ — -ii-i-. i i yj.o) -r \ — io,o4: i U-zuji-iC -f- ^u.zoo i u.uzuj^iic^ 

(-396.2 ± 0.4) + (0.90 ± 0.17)Ai + (0.170 ± 0.017)(At)^ 

— 139 ± 33 


S34, 15.5 

9009 07 


(302.3 ± 0.3) + (9.58 ± 0.07)At 


S35, 13.3 
2000.41 


(540.5 ± 0.1) + (1.70 ± 0.02)At 
(-437.8 ± 0.1) + (3.16 ± 0.03)At 


OOO, ID. 4: 

2004.56 


('97fif^-i-n9^-L(' 1 i'^-l-nifi'iA + 

\Z 1 O.O ZIZ U.Zj -f- l^ — 1,10 Hz U.lDJZAr 

(246.4 ± 0.3) + (-0.71 ± 0.17)At 


1 R 1 
OO/, ID.l 

2005.47 


yOOl..^ ZlZ U.^) -j- ^ — D.ID ZL. U.IO^/AC 

(390.2 ± 0.3) + (10.17 ± 0.13)At 


2006.94 


— 1/y.D ZlZ XJ.Xj) -\- y — OU.OD ZlZ U.DlJiiC -|- ^^ft.yZO ZIZ U.^toO jyLlZ) 

(55.2 ± 0.7) + (-10.91 ± 0.98)At + (-1.951 ± 0.754)(At)2 

]^g5 _|_ 


S39, 16.8 

^XJKJO- OKJ 


(-102.2 ± 0.7) + (-11.80 ± 0.16)At + (1.385 ± 0.132)(Af)^ 

/'ORQ Q -1- 1 91 _l_ f'^'i «1 -L 9Q1 At -U f 9 "^Sfi -1- 9'^^VA+I^ 


S40, 17.2 


(144.0 ± 1.5) + (3.96 it 0.48)At + (-3.772 ± 0.512)(At)2 

/oq 4 _|_ 9 Q\ _|_ /I 71 _|_ QS^A^^ -1- ( — 9 fi91 + 1 09fiVAf"\-^ 


S41, 17.5 


(-221.0 ± 0.6) + (-0.58 ± 0.37)At 

f 9QQ 0-1-0 fi'^ 4-/' 1QQ-I-0 "^81 Af 
^_ — zyy.u ^ u.Dj — i.yy ^ u.oo^iit 


S42, 17.5 
9004 QS 


(-160.1 ± 0.7) + (-6.18 ± 0.43)At 

T — "^^4 -1- 1 1 00-1-0 fi^1A + 


S43, 17.5 


(-493.2 ± 0.3) + (5.39 ± 0.12)At 

l"— 1'^4 4-l-0 4"14-f'7 7S-l-0 1R~lA^ 


S44, 17.5 
900R '^9 


(-92.2 ± 0.5) + (-9.48 ± 0.43)Ai 

('_94fi 0-1-1 0^-1- f' — 10 '^7 -1- 7fi^Af 


S45, 15.7 
onn'^ 47 


(193.3 ± 0.2) + (-6.61 ± 0.08)Ai 

( — ^-{^ -1- '^^ -1- ( — A Ofi-I- 1'^'lA/ 

^ i-O .\i ZIZ \).Oj 1^ '±.'i.)\j ZIZ U.-LO^i-iC/ 


S46, 15.7 

9001 4fi 


(246.1 ± 0.4) + (0.57 ± 0.12)Ai 
(—^7A -\- V\ A\ ( ^ '^7 -1- lO^lA^ 


S47, 16.3 

900R '^9 


(383.6 ± 0.8) + (-3.82 ± 0.70)Ai 

(OA^ 9 -1- 4^ -1- 0'^ -1- '^R^Af 


S48, 16.6 


(438.5 ± 0.5) + (-0.33 ± 0.14)At (-0.442 ± 0.098)(At)^ 

{4.70 1 -1- -1- n 9 9fi -i- 14"! A/ -\- ( — "^A"^ -\- OQSVA/I^ 


S49, 17.5 

900'=i fi^ 


(585.6 ± 0.6) + (15.61 ± 0.34)At 

/CI fi -1- 71 -1- n "^0 -1- '^4'1 A/ 


S50, 17.2 
200.^ 47 


(-504.7 ± 0.3) -1- (-2.84 ± 0.13)At 

/'_'^9« S -1- ^1 -1- fQ f^l -1-0 1R"\Af 


S51, 17.4 
200.^ 47 


(-473.5 ± 0.3) -1- (7.86 ± 0.14)At 

/'_9qQ 4 -1- ^1 -1- fS 01 -1-0 1 Af 
y — ^ifij.'* ZIZ ^ yo. yji. ZIZ \j,±oji^t 


S52, 17.1 
2006.94 


(200.8 -L 0.0) + ('2.68 -f- n.5fnA/ 


S53, 17.2 

9007 4R 


(323.7 ± 0.7) -1- (12.15 ± 1.04)At 

/ri^ 9 -1- 71 -U ('Q 04 -1- 1 0^1 A + 


S54, 17.5 

900fi i^iQ 


(135.4 ± 1.2) -1- (-1.05 ± 0.90)At 

/ fio 1 -1- Q"^ 4- (' 9R QO -1- R7^A + 


S55, 17.5 

900fi 7S 


(93.8 ± 2.9) (-23.72 ± 4.31)Af 

Qf^ R-l-'^ 01 -l-/'99 74-1-4 ISlAf 


S56, 17.0 

9007 SI 


(143.5 ± 1.2) + (-15.96 ± 2.81)At 
{ 1 '^4 fi-i-n^^l-i-i'nifi-l-i ifilA^ 


S57, 17.6 

9007 Aft 


(393.6 ± 1.1) + (-9.99 ± 1.43)At 

{ 1 47 A -\- f] 9i\ -L ( 4n'^-l-1 071 A f 
y — 111 .1 ^ u . 1 ^ I — ^ . u^) ^ 1 . u 1 1 z-i t 


S58, 17.4 
9nn'^ 47 


(-339.6 ± 0.4) + (5.98 ± 0.11)Ai + (0.189 ± 0.069)(At)^ 
{ f^fiQ 7-i-n 1'^1A+_L/'o i^s4 -I- osi^1 A +1^ 


S59, 17.2 
2007.89 


(-222.8 ± 0.4) + (5.87 ± 0.89)At 
(240.7 ± 0.8) + (-3.83 ± 1.71)At 


S60, 16.3 
2007.82 


(-277.6 ± 0.9) + (3.38 ± 2.19)At 
(168.0 ± 0.5) + (-19.17 ± 1.02)At 


S61, 17.9 
2007.90 


(-205.9 ± 1.4) + (-18.08 ± 2.17)At 
(-45.0 ± 1.4) + (-20.35 ± 2.12)Af 


S62, 17.8 
2007.90 


(-59.2 ± 1.3) + (-10.31 ± 2.34)At 
(67.8 ± 1.7) + (2.28 ± 3.47)At 
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t„ [vr] for (n. S) 

^{) [y^'\ ''z 


CK [mas] — 

S [mas] — 
u~ Lkiii/^J = 


S63, 17.5 
2007.93 


(196.0 ± 2.1) + (-12.04 ± 4.59)At 
(-134.7 ± 2.8) + (-9.44 ± 5.89) At 


Qfl/I 17 1^ 
iDCij 1 i .O 

2005.51 


f lOQ-Li 1^_L(' ip;/i9-i-n/iR^A + 

— IZ.O Hz ^-i-) ~r I, — ID.'iZ ZIZ \J.^K))i-l.Z 

(238.5 it 1.3) + (8.57 ± 0.57)At 


S65, 13.7 
onnn a i 


(-777.0 ± 0.1) + (2.17 ± 0.02)At 

( 9fiQ >=; -1- 1 -1- l" 1 4Q -1- n 09"^ A + 


S66, 14.8 
onnn 4i 

2004.63 


(-47.7 ± 0.2) + (12.90 ± 0.10)At + (0.015 ± 0.010)(At)^ 
( 1 nofi R -i- n 9 \ -L f 1 "iQ -4- n os"! _l /"n o47 -I- n ons^/' A+A^ 
12 ± 22 


QRT 19 1 
OD / , IZ, 1 

2000.41 
zuuo.yo 


41:01.41: i u.zj -|- ^ — lo.oy i u.u/ jiir -|- ^ — u.uod i u.uuo^^iic^ 
(872.4 ± 0.2) + (1.98 ± 0.08)At + (-0.058 ± 0.009)(At)^ 

1 =b 33 


S68, 12.9 
2000.41 


(275.8 ± 0.4) + (4.90 ± 0.07)Ai 
(764.9 ± 0.4) + (2.83 ± 0.08)At 


ooy, ±D.o 
2005.47 


(762.7 ± 0.5) + (0.69 ± 0.25)At 


o / u, iD.y 
2005.47 


( '^/tsi-l-ni'i-Lf' '?^'?-l-nni^^A-f 

— O-io. 1 i U.IJ -|~ — O.OO i \) AJO ) LXT 

(714.6 ± 0.2) + (-4.09 ± 0.08)At 


QVI 1 fi 1 
O ( 1 , lO. 1 

2005.47 
9nn7 AA 


(-762.1 ± 0.4) + (14.91 ± 0.09)Ai + (0.324 ± 0.068)(At)^ 


S72, 14.3 
onnn A^ 


(-668.1 ± 0.2) + (8.88 ± 0.03)At 


S73, 16.1 

9004 0*^ 


(-310.2 it 0.1) + (-9.93 ± 0.05)At 


S74. 16. n 
2005 47 


(-98.1 + n.2) + f-n,78 + n.09)A/ 

i' Nf; 1 1 1 II 9 ^ ^ J. fiJ 1 II (Kii \/ 

— ClUl.-i _L U.ZJ — |-i.U'd _L y> - \ — i L 


S75, 17.1 


(-154.5 it 0.2) + (6.03 it 0.12)At 
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